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Abstract 

This paper studies the iterates of the third order Lyness' recurrence X/.+3 = (a + 
Xk+i + Xk+2)/xk, with positive initial conditions, being a also a positive parameter. It 
is known that for a = 1 all the sequences generated by this recurrence are 8-periodic. 
We prove that for each o ^ 1 there are infinitely many initial conditions giving rise to 
periodic sequences which have almost all the even periods and that for a full measure 
set of initial conditions the sequences generated by the recurrence are dense in either 
one or two disjoint bounded intervals of R. Finally we show that the set of initial 
conditions giving rise to periodic sequences of odd period is contained in a codimension 
one algebraic variety (so it has zero measure) and that for an open set of values of a it 
also contains all the odd numbers, except finitely many of them. 
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1 Introduction and main results 



1.1 The third order Lyness' difference equation 

The excellent unpublished paper of Zeeman [12] about the celebrated Lyness' second order 
difference equation 

a + x n +i 

Xn+2 = , with a > , x\ > , > 0, (1) 

gives the key points for understanding the behaviour of the sequences generated by ([!]). In 
this reference it is proved that the map induced by (H|), 

f( x , y ) = L,^±y\ ( 2 ) 

defined on {(x,y) £M? : x > , y > 0} leaves invariant the level curves of the first integral 
V(x, y) = (x + 1) (y + 1) (a + x + y) /(xy) and, which is more important, that on each set 
{V(x,y) = h}, the map / is conjugated to a rotation of the circle with rotation number 
p a (h). By using this result, Zeeman explains the behavior of all the sequences generated by 
([1]), modulus a conjecture, the monotonous dependence of p a (h) with respect to h once a / 1 
is fixed. Recall that when a = 1, except for the fixed point, all the sequences generated by 
(PQ) are 5— periodic. This conjecture has been proved to be true in [2]. The study of the 
periods that can appear in the Lyness equation, as well as the study of the rotation number 
has also been done in pQ. 

This paper studies a similar problem to the one considered by Zeeman but in dimension 
three, and proves that in this case the dynamics are also described by rotations. The fact 
that we are in a higher dimension makes the problem more difficult. 

Concretely, we consider the third order Lyness' recurrence, 

a + x n+2 + x n+ i 

Xn+3 = , (3) 

for a > and positive initial conditions x\,X2 and £3, i.e. such that (xi, X2, x%) £ + : = 
{(x, y,z)€]R 3 : x>0,y>0,z> 0}. This recurrence is also known as Todd's recurrence, 
see [3 [10]. Recall that if some initial condition is such that (x±,X2, £3) = (xi+ P , X2+ p , xs +p ), 
and p is the minimal positive number satisfying this property it is said that this initial 
condition gives rise to a p— periodic sequence. It is well known that when a = 1 for any 
positive initial condition it holds that all the initial conditions in + of ([3]) are 8, 2 or 
1-periodic. We are interested to understand which is the situation when a / 1. Our main 
result is: 
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Theorem 1. Consider the third order Lyness' recurrence ([3jJ for a > and positive initial 
conditions x\,X2 and x%. 

(i) I/o^l there is a computable value qo(a) £ N such that for any q > qo(a) there exist 
continua of initial conditions giving rise to 2q-periodic sequences. 

(ii) The set of even periods arising when a £ (0, oo) contains all the even numbers except 
possibly 4, 6, 10, 12, 16, 18, 24, 28 and 40. 

(Hi) The set of initial conditions giving rise to odd periods is contained in an algebraic 
codimension one subset ofO + . Moreover, there is an open setlA C (0,1) U (1, oo) of 
values of a for which the set of the odd periods contains all the odd numbers except 
possibly finitely many of them. 

(iv) If a 1 then there exist a dense set of initial conditions in + such that the sequence 
generated by |5]) is dense in either one or two disjoint bounded intervals of R. 

The above theorem makes one to wonder the following natural questions: Are there some 
a > and some initial condition in + such that the recurrence at this initial condition 
is periodic of period 4, 6, 10, 12, 16, 18, 24, 28 or 40? Fixed a > 0, which are exactly all the 
even periods of the recurrence? And the odd periods? Is the set U introduced in Theorem 
□ (in) U = (0, 1) U (1, oo)? 

We want to remark that when the recurrence ([3]) is considered with initial conditions in 
the whole M 3 , the periods that can appear can be different. For instance in [3] it is proved 
that for some values of a there are initial (non positive) conditions giving rise to periodic 
sequences with periods 2,3,4,5,6,7 and Ap for any p > 3. 

The paper is organized as follows. In Section 11.21 we state our results on the discrete 
dynamical system generated by F thus obtaining the proof of Theorem [TJ All the results 
stated in Section 11.21 are proved in the following sections, moving some large proofs of 
the technical results to specific subsections and the appendices in order to improve the 
readability of the paper. 

1.2 Study from a dynamical systems viewpoint. 

As usual we reduce the study of the recurrence ^ to the study of the discrete dynamical 
system generated by the map F(x,y,z) = (y,z,(a + y + z)/x) defined in + . Note that 
this map is a diffeomorphism from + to + . A complete description of the discrete system 
generated by F gives a complete answer to the questions posed in Section II. 1\ and in 
particular a proof of Theorem Q] (see the end of this section). Our analysis of this dynamical 
system is done in two steps: 
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1. We will see that the phase space of F is foliated by invariant curves (sometimes 
degenerated to isolated points) which are given by the level curves of two functionally 
independent first integrals. The first step is to characterize the topology of this level sets, 
which turn to be diffeomorphic to circles (when they are not isolated points). 

2. The second step is to study the dynamics of F restricted to these invariant sets. As 
we will see, one of our main tools, at this stage, will be the study of an ordinary differential 
equation associated to the discrete dynamical system generated by F. This is an approach 
different to the ones in pQ, [12] (and even to the one in [2] although our starting point is 
the same of this last reference). Our approach turns out to be also effective for studying 
other difference equations, see [5]. 

Nevertheless there are some problems, named there as Questions 1 and 2, that have 
resisted our analysis. We remark that an answer to them would also allow to clarify the 
answers to the questions about ([3]) stated in Section 11.11 

Fixed a > 0, consider the diffeomorphism 



defined in 0+ := {(x, y, z) G M? : x > , y > , z > 0}. 

We begin by introducing some sets in + which are invariant under the action of F, in 
terms of the level surfaces of the well-known ([4, 6,[8l[TT]) couple of functionally independent 
first integrals of F, given by: 



Let L k = {(x,y,z) G + : V\(x,y, z) = k] and M h = {(x,y,z) G + : V 2 (x,y,z) = h} be 
the level surfaces of V\ and V 2 respectively. 

The orbits of F lie in 1^^ = n Mh for k > k c and h > h c , where k c and h c denote 
the values attached at the global minima in + of V\ and V2 respectively. For a given 
fixed h > h c , there exists k\ = ki(h),k 2 = k 2 {h) satisfying k c < k\ < k 2 and such that 
Ik,h 7^ only when k G [ki,k 2 ]. See Theorem [2] below, or Proposition [TT1 in Section [2] for 
more details about the topology of 1^. We will use the notation A = B to denote that the 
two manifolds A and B are diffeomorphic. 

Now, we introduce two interesting invariant sets that will play a very important role. 
The first one is 




(4) 



Vi(x,y,z) 



(x + 1) (y + 1) (z + 1) (a + x + y + z) 



V 2 (x,y,z) 



xyz 

(1 + y + z)(l + x + y)(a + x + y + z + xz) 



xyz 



C := {(x, (x + a)/{x - l),x) G M 3 such that x > 1} C O 
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It is easy to see that the set £ is a curve filled by two-periodic points of F and that it 
contains the unique fixed point in + : (x c , x c , x c ), where x c = 1 + \f\ + a. The second one 
is 

Q := {(x, y, z) G + such that G(x, y, z) = 0}, 

where 

G(x, y, z) = -y 3 - (x + z + a + l)y 2 - (x + z + a)y + xz(x + l){z + 1). (5) 

The set CuG is formed by the points in + where the gradients of Vi(x, y, z) and V^x, y, z) 
are parallel. 

In particular, it is not difficult to check that 

G(F(x, y, z)) = - a + y + Z G (x, y, z). (6) 
x z 

Note that this relation implies that Q is invariant by F and that F maps the zone {G > 0} 
into the zone {G < 0} and viceversa. Furthermore it implies that the dynamics of F 2 on 
the zone {G > 0} and on {G < 0} are conjugated, being the map F itself the conjugation. 
Figure 1 gives an example of the more generic position of L^, M^, Q and C. 




Figure 1: For a = 3, the level surfaces L31 and M42.5 which are diffeomorphic to spheres, 
the invariant surface Q and the line of 2-period points C. 
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Our main result about the dynamics of F, which proved in Section [3l is: 



Theorem 2. For each fixed h > h c . The following statements hold 

(i) Mh = Ufc e [fc li fc 2 ilfc ! / l , where the values ki, are given in Proposition\lll Moreover, for 
each k G (ki, A^), Ik,h splits into two disjoint connected components, that is: I^^h = ^kh^-^kh 
where I+ h := I kjh n {G > 0} = S 1 and I k h := I kjh n {G < 0} = S\ F(/±J = I* h , 
F 2 {Ikh) = ^kh an d the restriction of F 2 on each of these sets is conjugated to a rotation of 
the circle with rotation number p F 2(k,h). 

(ii) The set Q is invariant by F and the restriction of F to each set Q n {V\(x,y, z) = 

P 2 (k h(k)) 

k} , k > k c is conjugated to a rotation of the circle with rotation number pF(k) = — — 5 , 

where h(k) is a suitable known function. 

The proof of the above result is given in two steps. The first one is the study of the 
topology of the invariant sets I,, (this is done in Section [2]) . The second step is to prove 
that over these invariant sets F is conjugated to a rotation (this is done in Section [3]). The 
proof relies on some results that relate the rotation numbers associated to the invariant sets 
l£ h of F with the properties of a flow constructed from F which has the same invariant 
sets. 

Next results give some rotation numbers and periods appearing in the dynamical system 
generated by F. Their proofs use the regularity of the rotation numbers varying h, k and a. 
This regularity is studied in Section [H 

Theorem 3. For a > define 

1 / (l-a)vTTa 

p a := — arccos 



2vr \ v 2(l + ^TT^)(a + l + v / rT^) 

Then for each a / 1 there are circles of initial conditions in {G > 0} \ C and in {G < 0} \ C 
such that such that F 2 restricted to them is conjugated to a rotation with rotation number 
taking any value in (\, p a ) , if a > 1, and any value in [p a , \) if < a < 1. 

In Section [5] we give a constructive algorithmic approach to the problem of determining 
which are all the denominators of irreducible fractions which belong to a given interval, see 
Theorem [25] and Corollary 1261 In particular, by using these results and the above theorem, 
we prove: 

Corollary 4. (i) For any a^l there exists a computable value qo(a) £ N such that for any 
q > qo(a) there exist a continua of initial conditions giving rise to 2q-periodic orbits for F. 
(ii) Set 

'vr - 2 arcsin(l/8) I s 



£ rot 



4vr 3 
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Then, for each number p in I rot there exists some a > and a circle of initial conditions 
such that such that F 2 restricted to it is conjugated to a rotation with rotation number p. 
In particular, for all the irreducible rational numbers p/q £ I ro t, there exist periodic orbits 



(Hi) The set of even periods arising from the family {F(x,y,z) = (y, z,a + y + z/x) : 
a > 0} contains all the even periods except possibly 4, 6, 10, 12, 16, 18, 24, 28, and 40. 

The knowledge that we have of the odd periodic orbits of F is not so detailed as our 
knowledge of the even periods. We collect all our results in the following proposition: 

Proposition 5. (i) All the initial conditions giving rise to odd periods of F in + are 
contained in Q. 

(ii) There is an open set IA C (0, 1) U (1, oo) such that for each a £U the map F over Q 
has all the periods except possibly a finite number of them. 



For each p E J ro t there exists some a > and a circle of initial conditions contained in Q 
such that the map F restricted to them is conjugated to a rotation with rotation number p. 
Therefore, for all the irreducible rational numbers p/q £ J TO t, there exists periodic orbits of 
F of period q. 

All the above results and our numeric simulations of the functions pp2(k, k) and pp{k), 
detailed in Section [6l make as to propose the following questions. Note that the first one is 
similar to Zeeman's Conjecture. 

Question 1. Is it true that the function p F 2(k,h) varies monotonically when either k or 
h vary? 

If the answer is affirmative then all the rotation numbers, as well as the set of even 
periods given in Theorem [1] and Corollary [H are the only possible ones on + \ Q. 
Question 2. Is it true that for each a / 1, the rotation number pp is not identically 
constant on Q? Which is limit of pp when the initial conditions go to infinity over Q? 

For instance, the proof of Proposition [5] (ii) follows from the fact that, for a neigh- 
bourhood of values of a = ,„ 3 4 / c ° s ^^V v > , the rotation number on Q is not constant. Un- 

(2cos(27r/7)— l) z ' 

fortunately, we have not been able to obtain a general proof of this fact. Our numerical 
simulations for a = 3 and a = 7/9 (see Tables 1 and 3 in Section [6]) also show the same 
situation. If the answers to the above questions were affirmative we would obtain that in 
Proposition [5j U = (0,1) U (l,oo). 

The computation of the limit of pp over Q at infinity would give us useful quantitative 
information about which would be these odd periods of F. 



ofF* 



of period q. 



(Hi) Set 
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Remark 6. Note that when a / 1, each map F has infinitely many different periods and 
has sensible dependence with respect to the initial conditions. This last fact is because two 
close initial condition belong to two close sets, both diffeomorphic to circles, but over each 
one of them the rotation number is slightly different. 

Proof of Theorem Ql Parts (i) and (ii) are a direct consequence of Corollary [U Part (iii) 
follows from Proposition 

To prove (iv) observe that all initial conditions in + \ C give rise to rotations for F 2 
(respectively F). Moreover for most of these conditions, in the sense of Lebesgue measure, 
the associated rotation numbers are irrational. Therefore the orbits through these initial 
conditions are dense in a subset of Q (resp. + \ Q) which is diffeomorphic to S 1 (resp. 
the disjoint union of two S 1 ). The projection into the x-axis of the orbit of F coincides 
with the sequence generated by ([3]). This projection is formed by one or two disjoint closed 
intervals. Both situations are possible depending if the initial conditions are near the two 
periodic orbit or near Q. Hence the theorem follows. g 

2 Topology of the invariant sets of F. 
2.1 The results 

This section is devoted to prove the following weaker version of Theorem [2J Note that the 
difference between both results is that in this second one the dynamics of F or F 2 on each 
of the invariant circles is not yet described. The description of these dynamics is the goal 
of next sections. We use the following notations: A rti B means that A has a transversal 
intersection with B, and AU B means the union of A and B and that both sets are disjoint. 
Recall also that we say A = B when A and B are two diffeomorphic varieties. 

Theorem 7. (Topology of the invariant sets) Fix h > h c . Then 

(i) Mh = Ufc 6 [fc 1 fc 2 ]/fc,/i) where the values ki,k 2 are given in Proposition For each 
k G (^1,^2), Ik,h = ith ^ ^kh an( ^ eac ^ one °f these sets is diffeomorphic to a circle. 
Moreover F(/±J = % and F 2 (I± h ) = J± fc . 

(ii) The set Q is foliated by the fix point of F and the sets Q n {V\ = k} , k > k c , which 
are invariant by F and diffeomorphic to circles. 

The proof of the above theorem is done at the end of this section. To prove it, we first 
study the level sets of V\ and V2 in + and afterwards their relative position. 

Let L k = {(x,y,z) £ + : V\{x,y,z) = k} and M h = {(x,y,z) £ + : V 2 (x,y,z) = h} 
be the level surfaces of V\ and V2 respectively. It is well known that V\ has a global minimum 
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at (a; c , x c , z c ), where x c = 1 + \/l + a. We set 



(2 + VTTa) 3 (a + 3 + 3 \/TTa) 



(1 + \/rT^) 3 



Thus Lfc is not empty for k > k c , and Lk c = (x c ,x c ,x c ). 
Similarly, V2 also has a minimum at (x c ,x c , x c ). We set 




(l + yrr^) 



Then is not empty for h > h c , and M^ c = (x c ,x c ,x c ). 

Proposition [8] ( proved in Section [2.2p . states that except at the fix point all the level 
surfaces in + of V\ and V2 are diffeomorphic to spheres. Note that this result proves in 
particular that all the orbits of F starting at + lay in compact sets. 

Proposition 8. (General properties of Lp, and M^) 

(a) For k > k c , Lk is diffeomorphic to 8 2 . 

(b) For h > h c , is diffeomorphic to § 2 . 

Theorem [7] follows from the knowledge of the relative positions of the level surfaces Lk 
and M/j. First we describe the set T where and are not transversal and give the 
relative position of T and Lk- 

Lemma 9. (Locus of non-transversality of Lk and M^) Let T be the subset of + 
where Wi and W2 are linearly dependent, i.e T := {VV1IIVV2} H + . Then T = C U Q. 

Proof. Some computations show that 



(Vi) x (V 2 ) x 

(Vl)y (V 2 )y 



(z + 1) (1 + x + y) (a + z + y — xy) [ay 2 + ay — xz 2 — x 2 z+ 



y 2 + yz + y 3 + xy — xz + y 2 z — x 2 z 2 + xy 2 ) / (x 3 y 3 z 2 ) 



{Vi) x (V 2 ) x 
[V X ) Z (V 2 ) z 



(y + 1) (x — z) (a + x + y + z + xz) [ay 2 + ay — xz 2 — x 2 z+ 



9 q o 9 9 9\//^9^\ 

y + yz + y + xy — xz + y z — x « + xy ) /{x y z ) 



and 



(Yi) z (v 2 ) z 



(x + 1) (1 + y + z) (a + x + y — yz) (ay 2 + ay — xz 2 — x 2 z+ 

9 S 9 9 9 9\ // S S 9\ 

y + yz + y + xy — xz + y z — x z + xy J /(x y z ). 
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The solutions in + of the above three functions equated to zero satisfy either 

ay 2 + ay — xz 2 — x 2 z + y 2 + yz + y 3 + xy — xz + y 2 z — x 2 z 2 + xy 2 = 0, 

which is precisely Q or {(x, y,z) : a + y — xy + z = 0, x — z = 0, a + x + y — yz = 0} which 
coincides with C = {(x,y,z) : y = (x + a)/{x — 1), z = x}, as we wanted to prove. g 

The topology of T is given by the next result, which is proved in Appendix [Cl 
Proposition 10. (Topology of the nontransversality locus) Fix k > k c . Then 

(i) C n Lk consists of two points which are a 2-periodic orbit of F . 

(ii) G (ti L k . 

(m) gni^s 1 . 

To describe the relative positions of the level surfaces and Mh we keep with 
h > h c fixed and consider for all k > k c , obtaining: 

Proposition 11. (Relative positions of and M^) For a given fixed h > h c , there 
exists k\ := k\{Ji),ki := ki(K) satisfying k c < k\ < k% and such that the following statements 
hold: 

(a) If k € [fcofci) then I k ,h = 0- 

(b) If k = k\ then either It x ,h = C D Mh (which are two points describing a 2-periodic 
orbit), or I klA = GnM h ^S l . 

(c) Ifk G (h,k 2 ) then I kyh = S X U S 1 . More precisely, h,h^{G > 0} = S 1 and I k)h Pi{G < 
0} S S 1 . 

(d) For k = k 2 then either I hih = CnM h if I kl>h = Q n M h , or I k2>h = Q D M/, = S 1 , i/ 

i kl , h = cnM h . 

(e) Ifk>k 2 then I h)k = 0. 

The proof of Proposition [TT] is given in Subsection 12.31 Now we can prove Theorem [3 

Proof of Theorem^ The result follows directly from Proposition [11] and the above explained 
consequences of expression j6|). 1 
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2.2 Proof of proposition [8] 

To study the surfaces solving Vi(x,y, z) = k, we get that they can be written as the 
union of the graphs of the two functions Z- and z + , given by: 



a(x,y;a,k) ± y/A(x, y\ a, k) 

z±{x, y; a, k) = — r- , (7) 

0{x,y;a,k) 

defined in {(x,y) G R 2 : A(x,y; a, k) > 0}, where 

a(x, y; a, fe) = —a — 1 — (a + 2) x — (a + 2) y — x 2 — (a — k + 3) xy — y 2 — x 2 y — xy 2 , 

(3(x, y;a,k) = 2(1 + x + y + xy), 

and 

A(x, y; a, k) = (a - l) 2 + 2a (a - 1) x + 2a (a - 1) y + (2a - 2 + a 2 ) x 2 

+ (-2k - 2ka + 4a 2 - 2) xy + (2a - 2 + a 2 ) y 2 + 2ax 3 + (-4/c - 2ka - 2 + 6a + 2a 2 ) x 2 y 

+ (-4A; - 2/ca - 2 + 6a + 2a 2 ) xy 2 + 2ay 3 + x 4 + (4a - 2fc + 2) x 3 y 

+ (-2/ca - 6k + a 2 + 3 + 6a + k 2 ) x 2 y 2 + (4a -2k + 2) xy 3 + y 4 + 2x 4 y 

+ (-2A + 4 + 2a) x 3 y 2 + (-2A; + 4 + 2a) x 2 y 3 + 2xy 4 + x 4 y 2 + 2x 3 y 3 + x 2 y 4 . 

Observe that 

(a 2 - A)(x, y; a,k)=A(y + I) 2 (x + l) 2 (x + y + a) > 0, (8) 

for (x,y) G Q + := {(x,y) : x > 0, y > 0}. Hence z±(x,y; a, k) ^ on Q + . This means that 
either z±(x, y; a, k) £ + for all (x, y) G Q + , or z±(x, y; a, k) G 0~ := {(x, y,z):x>0,y> 
0, z < 0} for all (x,y) G Q + . In particular, each connected component of with x > 
and y > is completely contained either in + or in 0~ for all k > k c . 

On the other hand notice that each level surface has an equator described by 

z±(x, y; a, k)\ A ^ x y . a ^ =0 . 

A description of the planar algebraic curves := {(x, y) G M 2 : x > 0, y > 0, A(x, y; a, k) = 
0} is given in the next lemma, which will be proved in Appendix O See Figure 2 for more 
details. It is the key result to prove Proposition [8] (a). 

Lemma 12. For k > k c the planar algebraic curve := {(x,y) G M? : x > 0,y > 
0, A(x, y; a, k) = 0} consists of 

(a) If k > k c : two concentric ovals 7^ and £fc surrounding the point (x c ,x c ). Furthermore 
7^ shrinks to (x c ,x c ) when k — > k c , and if a < 1 then the oval Cfc has a contact with 
the axis {x = 0} and {y = 0} at (0, 1 — a) and (1 — a, 0) respectively. 
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(b) If k = k c : one oval Q kc o,nd the point (x c ,x c ) G Int(Cfc)- 



y 




X X 

Figure 2: The curve of Lemma [T2l 



Proof of Proposition^ (a). By Lemma fT2| for any k > k c , Q + is split in the regions A k , B k 
and C k , as is shown in Figure 2, defined in the following way: 

A k = Int(7 k ), 

B k = Int(C k )\{7kUlnt( 7k )}, 
C k = Q+\ {(k Ulnt(Ck)}. 

Now we will see that for any k > k c we have A(x, y; a,k) > for all (x, y) G U C k \ 
and A.(x,y;a,k) < for all G B k . This means that the surface {V\ = k}, defined by 

([7]) only exists for (x,y) G A^U C k - 

Indeed, we can write A(x c , x c ; a, k) = x 2 y 2 k 2 + pi(a)k +po(a), and on the other hand 
A(cc c , x c ; a,k) = for k = k\ and k = k c , such that k\ < k c . Hence, for k > k c we have 
A(x c , x c ; a, k) > 0, and this proves that A(x,y; a, k) > for all (x,y) G A k . On the other 
hand A(0, 0; a, fc) = (1 — a) 2 > 0, hence A(x, y; a, fc) > for all (x, y) G C k . 

Finally, it can be seen that the zeros of A(x, y; a,k) = on Q + are simple so that 
A(x, y; a, k) < for all (x,y) G B k . 

Now we observe that lim y _ +0 + a(x, y; a, k) = a(x, 0; a, k) = —x 2 — (a + 2) x — a — 1 < 0, 
and that lim a(x, y; a, k) = —oo, for all x > and fc. This means together with the above 

y— >+oo 

observation concerning equation ©, that {Vi = fc} C O - , for all (x,y) G Cfc. 

Observe that z±(x c ,x c ;a,k c ) = x c > 0, hence by continuity z±(x c , x c ; a, k) > for 
k > k c . But, as seen before, by equation (jSJ), each connected component of {V\ = k} 
with x > and y > is completely contained either in + or in 0~ for all k > k c . So 
2±(«,?/;o,A;) > for (x,y) G A k . 
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Therefore is given by ([7]) for (x, y) G hence is a topological sphere. To see that 
indeed it is diffeomorphic to a sphere, by using the implicit function Theorem, it suffices to 
prove that the function V\ has no critical points on Computing the partial derivatives 
of Vi we get: 

dVi (y + l)(z + l)(-x 2 + a + y + z) 

dx x 2 yz 

dVi (x + l)(z + l)(-y 2 + a + x + z) 

dy xy 2 z 
dVi _ (x + l)(y + l)(-z 2 + a + x + y) 

dz xyz 2 
Hence the critical points of V\ which lie on + have to satisfy x 2 = a + y + z , y 2 = 
a + x + z , z 2 = a + x + y which easily implies x = y = z = x c . So the only critical point of 
V\ on + is the fixed point. Hence part (a) of the proposition follows. g 

To prove Proposition [5] (b) we proceed in a similar way that in case (a). Solving 
V"2(x,y, z) = h, we get that the surface {V2 = h} can we written as the union of the 
graph of the two functions: 

z±(x,y;a,h) = ^ ± V^ Kg , ( 9 ) 

P{x,y,a,h) 

where 

a(x, y; a, h) = —ya — 2x 2 — 3x — a — y 2 x — yx 2 — 5yx — 3y + kxy — xa — 2y 2 — 1, 

0(x, y; a, h) = 2 + 2x 2 + Ax + 2y + 2yx, 

and A (a;, y;a,h) = 1 — 2a + 2x + 2y + h 2 x 2 y 2 — 6yxa + 2ya 2 x — Ay 2 ax — 2y 3 ax — Ay 2 ax 2 — 
4yax 2 — 2x 3 ya — 6hx 2 y — 6kxy 2 — I0hx 2 y 2 — 2hx 2 y 3 — 2hx 3 y 2 — 4hx 3 y — 4kxy 3 + 4yx 2 + 
a 2 — Axa + Ayx — 4ya — 2kxy + x 2 + y 2 + Ay 2 x + 2xa 2 — 2kxy 2 a — 2hx 2 ya — 2kxya + y 2 a 2 + 
2x 2 y 3 + 2xy 3 + 2ya 2 — 2x 2 a + 2x 3 y + x 2 a 2 — 2y 2 a + y 4 x 2 + y 2 x 4 + 2x 3 y 2 + 5x 2 y 2 + 2y 3 x 3 . 

By looking at the above coefficients it is easy to check that, if x > and y > 0, 
a(x, y; a, h) 2 — A > and hence z±(x, y; a, h) ^ on Q + := {(x, y) : x > 0, y > 0}. Thus, 
either y; a, /i) G + for all (x, y) G Q + , or y; a, /i) G O - := {(x, y, z) : x > 0, y > 
0, z < 0} for all (x, y) G Q + , that is each connected component of {V2 = h} with x > and 
y > is completely contained either in + or in 0~ for all h > /i c . 

On the other hand observe that each level surface has an equator given by the equa- 
tion z±(x,y;a,h)\ A r xy . ah \_ . The description of the planar algebraic curves, dropping the 
subindex a, Th := {(x, y) G M 2 : x > 0, y > 0, A(x, y; a, h) = 0} is again the key of the proof 
of Proposition [5] (b), see Figure 3. We will use the following lemma, proved in Appendix [Bl 
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X 

Figure 3: The curve Th of Lemma [T3"l 

Lemma 13. For h > h c the planar algebraic curve Th := {(x,y) G M 2 : x > 0, y > 
0,A(x,y;a,h) = 0} consists of 

(a) Two branches y = yi(x) , i = 1,2, such that yi(x) < y2(x), for x > 0. Furthermore 
these two branches satisfy lim yAx) = +oo and lim ydx) = + . 

(b) An oval^fh, contained between these to branches ifh > h c , and a single point ifh = h c . 

Proof of Proposition [3 (b). By using Lemma [T3l we have that for any h > h c , Q + splits in 
four regions A^, B^, Ch and D^, as is shown in Figure 4, defined in the following way: 

A h = {(x,y) G Q+ :y< yi(x)}, 
B h = {{x,y) eQ + :y> y 2 (x)}, 
Ch = 7/ l Ulnt(7 h ), 

D h = {(x,y) GQ+ <y<y 2 (x)}\{ 7 ^Ulnt(7h)}. 

It is easy to check that for any h > h c , A(x, y; a,h) > for all (x, y) £4UB/,U Ch', and 
A(x,y; a, h) < for all (x,y) £ Dh- This means that the surface {V2 = h}, is only defined 
by © on the region A h \JB h llC h . 

We now observe that lim y ^ + a(x, y; a, h) = a(x, 0;a,h) = — 1 — a — (0 + S)x — 2x 2 < 0, 
and lim a(x,y;a,h) = —00, for all x > and h. Hence {V2 = h} C 0~ , for all (x,y) G 

y— >+oo 

4UB/,. 

We observe that 2±(a; c , x c ; a, /i c ) = x c > 0, hence by continuity z±(x c , x c ; a, h) > for 
/i > /i c . But, as seen before each connected component of {V2 = h} with x > and y > 
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is completely contained either in + or in O for all h > h c . So z±(x,y;a,h) > for 
(x,y) e C h . 

Therefore Mh is given by ([9]) for (x,y) £ Ch, hence Mh is indeed a topological sphere. 
Finally, let us see that the surface Mh is a differentiable manifold for h > h c . It is enough 
to see that V2 has no critical points on Mh- The partial derivatives of V2 are: 



Qy 2 _ (y+z+l)( y -x 2 -x 2 z+ya+y 2 +yz+a+y+z) 
dx x 2 yz ' 

dV2 a+x+z+xza+xa+x 2 +'i xz—2y 2 —xy 2 z+za+z 2 +2x 2 z—2xy 2 —y 2 a—2y 2 z+2 xz 2 +x 2 z 2 —2y s 

dy xy 2 z ' 

QV 2 {x+y+l)(-z 2 -xz 2 +ya+xy+y 2 +a+x+y) 
dz xyz 2 ' 

The critical points on + have to satisfy: 
p(x, y, z) := —x 2 — x 2 z + ya + y 2 + yz + a + y-\-z = 

q(x, y,z) := a + x + z + xza + xa + x 2 + 3 xz — 2 y 2 — xy 2 z + za + z 2 + 2 x 2 z — 2 xy 2 

-y 2 a - 2 y 2 z + 2 xz 2 + x 2 z 2 - 2 y 3 = 
r(x, y, z) := — z 2 — xz 2 + ya A- xy + y 2 + a + x A- y = Q. 

Since p(x, y) — r(x, y) = (z — x)(l + x + y + z + xz) we get z = x, and substituting this 
equality in p(x, y, z) and q(x, y, z) we have the system: 

s(x, y) = —x 2 — x 3 + ya + y 2 + xy + a + y + x = 0, 

t(x, y) = (x + y + 1) (x 3 + 3 x 2 — yx 2 + xa + 2x — 2xy — 2y 2 + a — ya) = 0. 
Denoting 

. . Q Q Q Q 

u(x,y) = x +3x — yx +xa + 2x — 2xy — 2y +a — ya 

we get 

2s(x, y) + u(x, y) = — x 3 + x 2 — yx 2 + xa + 4x + 3a + ya + 2y = 

which let us to isolate y in terms of x : 

x 3 — x 2 — (4 + a)x — 3a 
V ~ -x 2 + a + 2 ' 

Then 

/ x 3 - x 2 - (a + 4)x - 3a\ h(x) 
s[x,- 



(10) 



(-x 2 + a + 2) ; (_ x 2 + a + 2) 2' 

where a(x) = — (a — l + x + x 2 ) (a + 2 x — x 2 ) (2a + xa — x 3 — 2 — 2x — 2 x 2 ) , and hence 
we have to consider three cases depending on the zeros of h{x). 

If a — l + x + x 2 = and x > 0, then x = x c and from (fTUI) . y = x c . Since z = x we get 
the fixed point. 
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Ifa + 2x — x 2 = and x > 0, then x = (— 1 + — 4o)/2 and substituting this value of 
x at (fTUj) we see that the corresponding y is negative, so we do not need to study this case. 

Now we have to consider the positive roots of g(x) := x 3 + 2 x 2 + (2 — a) x + 2(1 — a). 
We notice that when a < 1 there are not changes on the signs on the coefficients of g(x), 
and hence there are not positive roots of g(x) = 0. When a > 1, then there is a unique 
change of signs between the coefficients of g(x), and hence, by the Descartes rule, and since 
#(0) < and lim g(x) = +oo, we get exactly one positive real root, say x. We claim that 
for x = x, the corresponding value of y given in (|10p is negative. To prove this observe that 
g{Va + 2) = 4Va + 2 + 4a + 6 > 0, which implies that x < \Ja + 2. So, a + 2 - x 2 > 0, i. 
e., the denominator of (|10p is positive. By evaluating the numerator of (|10p at x, we get: 

x 3 - x 2 - (4 + a)x - 3a = x 3 - x 2 - (4 + a)x - 3a - g{x) = -(3x 2 + 6x + a + 2) < 0. 

Hence, there are no critical points of V2 in + different from (x c ,x c ,x c ), and the result 
follows. I 

2.3 Proof of Proposition [TT1 

To prove Proposition [TT] we need some technical results, stated below. 
Lemma 14. (a) I kcihc = (x c ,x c ,x c ), and I k>hc = %fork> k c . 
(b) Fix h > h c . Set k\ = minVi and ki = maxl^, then the following statements hold. 

M h M h 

(i) h,h / if and only if k G [fci, k 2 }- 
(ii) Ik- h C T = Q U C, for each i = 1,2. 

Proo/. (a) For k = k c L kc = M hc = (x c ,x c ,x c ), hence h c ,h c = (x c ,x c ,x c ). For k > k c , since 
(x c ,x c ,x c ) £ L k , I k>hc = 0. 

(b) Since Mu is compact then there exists k\ = min V\ and ko = max Vi . 

M h ~ M h 

Observe that ki 7^ k\ because V\{M^) is not constant, since otherwise 
VV 1 (M h )\\W 2 (M h ), but this only happens in M h n (£ U Q) / M ft (to prove this last 
inequality just consider that if h > h c , then D £ consists of two points which are not 
contained in G, hence Mh n (£ U Q) is disconnected). 

By definition of k\ and A:2 it is obvious that I k h = if A; ^ [k\, k 2 ]- 

Let a, 6 G M 3 be a points in Mh such that Vi(a) = fej, and Vi(6) = ^2 (observe that this 
points exist because the absolute extrema of Vi\{M h } are reached). 

Take now a continuous curve 7 : [0, 1] — * Mh such that 7(0) = a and 7(1) = b. The 
function 

a: [0,1] — > [fe^fca] 
f _^ Vi( 7 (t)) 
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is continuous, that is for all k 6 [k\,k2) there exists at least t k such that g{t k ) = k, and 
7(ifc) G Ife^. Hence the result follows. 

(ii) By the definition of k\ and k2, and using the theory of extrema with constraints, V\ 
reaches these values at points where the gradient vectors of V\ and V2 are parallel, hence in 
J 7 , as we wanted to prove. g 



Lemma 15. Fixed k > k c . ^2|{£ fc n£} = Cl anc ^ ^ 2 \{L k ng} = C2 w ^ iere °i an d c 2 are different 
constants. 

Proof. From Lemma [TOl (i), L k n £ = p\ U p 2) but since {^1,^2} is a 2-periodic orbit and 
V 2 is an invariant V^pi) = V2(F(p 2 )) = ^2(^2) = c\. 

Also V2 is constant over n Q. Indeed, by Lemma flUl (ii), L k n Q = S 1 , hence we 
can consider a (^-parameterization of L k n <7 given by 7(i). By definition ^^i(7(*)) = 
VVL(7(t))-V(t) = 0. But observe that on 0, Wi||W 2 , hence W 2 (7(*)) VC*) = Wi(7(t))- 
= 0, thus V2(7(i)) is constant. Hence V-^Lfc D (7) = C2. 

To prove that c\ 7^ C2, just observe that by a same argument than the one used in the 
proof of Lemma [TH c\ and C2 both must be the extrema of V^Lfc), and V2 is not constant 
over L k . | 



Corollary 16. (%) Either 

h u h = CnL k ^piU p 2 , and I k2>h = £ n L fc = S 1 , or 
hx,h = Q n L fc = S 1 , and J fc2)ft = £ n L fc = pi U p 2 , 
where pi,p2 is the two-periodic orbit, located in L k . 
(ii) If k G (fei, A; 2 ) t/ien iti M^. In particular I k h = U fi n it e *§ ■ 

Proof, (i) Lemma [HI (ii) ensures that I kl h G .T 7 = £ U £/. Lemma [15] prevents that there 
exist points a,b £ L kl n such that a € £ and b £ Q. Thus either we have 

• Case 1: I kljh C £ and hence = I klth nC, or 

• Case 2: J fcl)/l C £/, (hence J fel) ^ = J feljft n £). 

Observe that by Lemma [TU1 (i) Lfc a n £ = p\ U P2, but since, by Lemma [T5l V^pi) = 
V2 (p2 ) j we have that in the Case 1 , I kl j t = L kl n M\ t n £ = L kl n £ = pi U p2 • 

In the second case, we need to prove that Ik lt h = L kl D Pi <5 = n = S 1 . Observe 
that this is a consequence of the fact that by Lemma [T5l V2 is constant over L kl n £/. 

The same argument holds if instead of I kl h we consider I k2 ,h- But now observe that if 
we are in the first case, then I k2) h = L k Pi Q, because each point belongs only to one level 
set of V\. The same happens in the second case. 
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(ii) Statement (i) implies that the locus of non transversal intersections of the foliation of 
+ given by {Lk}{k>k c } with Mh are given only by I^hi i = 1; 2. On the other hand Lemma 
[T4l(i) ensures that Ik,h 7^ for k G {k\, k 2 \ therefore L k and Mh must intersect transversally 
in this case. Thus, see [HI page 30], Codim(Lfc C\Mh) = Codim(L/ c ) + Codim(M/ l ) = 1 + 1 = 2 
and Lfc Pi Mh is a submanifold of M 3 . This implies that I k h is a union of curves. But since 
both Lfc and Mh are compact, and each connected, compact 1 -dimensional manifold is 
diffeomorphic to S , see page 208], then I kh = US 1 . But these disjoint union of S 1 lie 
in a compact region (say Mh) and are defined by analytic equations, therefore it must be a 
finite union. g 

Next lemma shows that in (ii) of the above Corollary the finite union is exactly two § . 
Lemma 17. (Topology at the transversal intersections of L k and M/J For all 

ke (h,k 2 ), l kth ^S 1 u s 1 . 

Proof. From Corollary [TU] (ii) and Proposition [5] (b) we know that for all k £ (ki,k 2 ), 

h,h — UfiniteS 1 - 

Observe that Vi = {-^A;}{fee(*i,fe)}> induce a foliation of closed curves on Mh nesting the 
2-periodic points defined by Mh Pi C. These two periodic points are the only ones in the 
foliation of Mh induced by V2 = {-^jfc}{fee[fc 1 ,fc 2 ]}) an d since they are in the plane z = x, each 
of the closed curves of Vi must intersect the plane z = x. We want to prove there are only 
two of them. 

Consider now the restriction of Mh to the plane z = x, given by the equation V2(x, y, x) = 
h, which solutions are described by two functions x — » y±(x, h). 

Let vi(x) := V\(x, y+(x, h),x) be the restriction of V\ over the branch y = y+(x, h). We 
only need to prove that fixed k G (ki, k 2 ), the equation 

V 1 (x,y + (x,h),x) = k, (11) 

has only two solutions, that correspond to the two closed curves of the statement (observe 
that from expression for any closed invariant curve 71 such that 71 Pi {G < 0} 7^ 0, we 
have 71 n {G > 0} = 0). To see this we will prove that the singular points of v\(x) are 
located in (Q U C) n {z = x}, hence v\ is monotonic for those x such that (x, y+(x, h),x) E 
{G > 0} D {z = x} or {G < 0} D {z = x}, and therefore equation (jlip has two solutions. 
Indeed, in z = x, (Vi) x = (Vi) z for i = 1, 2, 



v[(x)= ^(V 1 ) x + (V 1 ) y ^± + (V 1 ) z 



{ z=x ,y=y + } 



2(V 1 ) x + (V 1 ) y dy 



dx 



{z=x,y=y+} 

(12) 
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and from V 2 (x, y, x) = h, we have 

dy+ f(V 2 ) x + (V 2 ); 



2(V 2 ) X 

{z=x,y=y+} \ ( V 2)y 



(13) 

{z=x,y=y+} 



dx \ (V 2 )y 

Using equations (fT2|) and (fT3|) we have v' 1 (x) = if and only if 

(Vl)x _ (Vl)y 

(V 2 ) X {V 2 )y 

Hence on the locus where the gradient vectors of V\ and V 2 are parallel. This set is (£ U 
Q) n {z = x}, as we wanted to prove. g 

Proof of Proposition ITTl Statements (a) and (e) are a direct consequence of Lemma [1~4"1 
(i), statements (b) and (d) from Corollary 1161 (i) and, finally, statement (c) follows from 
Lemma [TTJ g 

3 Dynamics of F. Proof of Theorem [21 

Next result relates, under some hypotheses, the dynamics of an ordinary differential equation 
and a discrete dynamical system that share an one dimensional invariant set. 

Theorem 18. Let f : U — » U be a C 1 map where U is an open connected set U C M n , 
satisfying the following assumptions: 

(Al) There exists a C 1 vector field X in U such that 

X(f{q)) = (Df) q X(q) for all q G U. (14) 

(A2) For a fixed p G U, the map f leaves invariant j p , where ^ p is the orbit of x = X(x) 
which passes trough p. In particular there exists r£l such that ip(r,p) = f(p), where 
(f(t,p) is the flow of X. 

Then, if j p = S , the restriction of f on -y p is conjugated to a rotation on the circle with 
rotation number t/T, where T is the period of ' 7 P . 

Proof. By substituting q = <p(t,p) in (fl~4"]) we obtain: 

X(f^(t,p))) = (Df)^ p) X( V (t,p)) for all t G R. 

Notice that the above equality precisely says that the function t — > f(f(t,p)) is also a 
solution of x = X(x). Since when t = it passes trough f(p), by the theorem of uniqueness 
of solutions we have that 

/Op)) = f(v{t,p)) for all ie/ p n i m . 
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Hence, since ip(r,p) = f(p), if q = p(t,p) we get that ip(t,ip(r,p)) = /(g) or, equivalently 
that 

c p(r,q) = f(q) for all <?G7 P . (15) 

Let us prove by using this relation that the map / : 7 P — * j p is conjugated to a rotation of 
the circle with rotation number p := r/T. 

Indeed we prove that the map h : S 1 — > 7 P given by /i(exp (■££)) = <-p{w^t,p\ is the 
desired conjugation. To see this it suffices to show that / o h = h o ry, where 77 is the 
rotation of angle 2ttt/T. The following chains of equalities give us the desired result. 

(/ h) (exp (it)) = f (u> (J^t, pX) =<f(r+ ^t, >> 

(hor T )(exp(it)) = h (exp (i(t + ^)) J = <p ( ^(t + ^r),p) = vf^ + ^^/ f 
where we have used (TT5lh 



Til r V 2vr v T " / r V 2?r 



Proof of Theorem By using Theorem [TJ only remains to prove that F or F 2 restricted to 
the invariant leaves given in this theorem are conjugated to rotations. This will be done by 
using Theorem [TSl To apply it we need a vector field X having the same invariant leaves 
that in Theorem [7] and satisfying (|14p . We start with the vector field X = Wi x VV2, 
where recall that V% and V2 are the invariants of F. We obtain that 



Xi(x,y,z) 
X 2 {x,y,z) 
X 3 (x,y,z) 



(x + 1)(1 + y + z)(yz — x — y — a)G(x, y, z)/(x y z ), 
(y + l)(z — x)(a + x + y + z + xz)G(x, y, z)/(x 3 y 2 z 3 ), 
(z + 1)(1 + x + y)(a + y + z - xy)G(x, y, z)/(x 3 y 3 z 2 ). 



Clearly it has V\ and V% as first integrals, but unfortunately it does not satisfy (JHJ). It 
is natural to try to remove the common factors of the components of the above vector 
field. We consider the new differential equation defined by the vector field X(x,y, z) := 
^^(VVi&y^) xVV 2 (x,y,z)): 



x = Xi(x,y,z) 
y = X 2 (x,y,z) 
z = X 3 (x,y,z) 



(x + 1)(1 + y + z)(yz - x - y - a)/ (yz), 
(y + l)(z — x)(a + x + y + z + xz)/(xz) (16) 
(z + l)(l+x + y)(a + y + z- xy)j (xy), 

A computation shows that it satisfies condition ([ED, i.e. X(F(q)) = (DF) q X(q) in + , 
and then also X(F 2 (q)) = (DF 2 ) q X(q) in O+. 

Since X also has Vi and V2 as a first integrals, each connected component of Ik t h will be 
an orbit of x = X(x). By Theorem the sets I kjh n {G > 0} = S 1 and I kjh n {G < 0} = S 1 , 
for k £ (k±,k2), are periodic orbits of X and invariant by F 2 . Since condition (|14|) is 
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satisfied for F 2 , Theorem [18] applies and F 2 is conjugated to a rotation of the circle. Hence, 
statement (i) follows. 

(ii) Now consider Q n {Vi = k}, which by Theorem[7]is also a periodic orbit of x = X(x) 
and is invariant by F. Since (|14p is also satisfied, by using again Theorem [18] we get that 
on Q n {V\ = k}, F is conjugated to a rotation of the circle, as we wanted to prove. g 

4 Properties of the rotation numbers of F 

The main result of this section proves the analyticity of the rotation number of F 2 in 
+ \ {£} and computes the limit of these rotation numbers when we tend in a certain way 
to the line of two periodic points C 

Proposition 19. (i) For each fixed a > and h > h c consider the values k\ := k\(a) 
and &2 := ^(a) given in Theorem I 111 Then the assignment (k,h,a) — ► p(k,h,a), where 
p(k,h,a) is the rotation number of F 2 restricted to Ik,h, is analytic for all (k,h,a) with 
a > 0, h > h c and k 6 (k\, Z^)- 

(ii) Fix a > and k > k c . Let pF{k) be the rotation number of F at each level curve 
{Vi = k} n Q, and let —1, cos(^) ± i sin(^) be the three eigenvalues of DF at the fix point of 
F. Then 

v /t \ 1 {a-l + VT+^L\ 

lim PF(k) = — = — arccos . 

k^k c r y ' 2ir 2vr \ 2a / 

(Hi) Fix a > and consider the surface Mh with h > h c a fixed value. The set Mh(~){G > 
0} is filled by closed curves given by jrj". and by a fix point of F 2 given by p^ = Mh n {V\ = 
k*} n {G > 0} € £ given by ph = (xh, (a + — l),Xh), x^ > 1. Let PF^{k) be the 

rotation number of F 2 on I~£ h and 1, cos (0ft ) ±zsin(0ft) the three eigenvalues of DF 2 at ph- 
Then 

k-^k* 2tt 2tt \ 2xft(a + XftJ / 

To prove the above proposition we need some preliminary results. 

Lemma 20. (i) The only singular points of the vector field $16\) in + are on C Moreover, 
except on C, the plane E = {z = x} is a global transversal section for its flow. 

(ii) Set a = ao > 0, h > h C) k G ^2), and po = (xo, yo, xo) G E \ {£ U Q}. Then there 
exists a neighborhood of po in E (namely Ej oc = E n B e (po) ), such that the points p G E^ oc 
depend analytically on a, k and h. 

Proof. A straightforward computation shows that X is orthogonal to E outside C and that 
C is filled by the only singular points of X in + . So, statement (i) follows. 
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It is important to notice that as we will see in the proof of Lemma [T71 all the periodic 
orbits of the vector field X must intersect S. 
To prove (ii) consider 

V : M 5 — ► M 2 

(x,y,a,h,k) — ► (Vi(x,y,x) - k,V 2 (x,y,x) - h), 

where the dependence of a is hidden in V\ and V%. On one hand V(xq, yo, ao, &o> ^o) = 0) 
and on the other hand 

(x + l)(x + y + l)(a + x - + y) G(x, y, x) 

X 5y3 

z=x 

which is nonvanishing in £ \ {£ U £7}. From the implicit function theorem, x and y are 
analytic functions of a, k and /i in a neighborhood of po = (xq, yo, xq), namely S; oc . | 

Lemma 21. Ze£ <p(t,p, a) be the flow of [T6)) . where we explicitly write the dependence with 
respect to a. Fix a > and take p = (x,y,x) £ Y>\{£{J Q}. 

(i) If T(jp,a) is the period of the periodic orbit of X passing through p, then T(p,a) is 
analytic at (p, a) . 

(ii) Ifr(p, a) is defined by the equation <p(t(p, a),p, a) = F 2 (p), then r(p, a) is also analytic 
at [p, a). 

Proof. Consider the system (p(t,p, a) — p = 0. Obviously ip(T(p, a),p, a) — p = 0, and 
^ (<p(t,p,a)-p) = X(p) = (Xi(p),0,-Xi(p)) + 0, 

because X\ (p) ^ in By applying the implicit function theorem to the first component 
of the above system we have that in a neighborhood of (p, a) in X\{£U£?} the period function 
T(p, a) is analytic. The proof of (ii) follows applying the same argument to equation 
<p(t,p,a) -F 2 (p) = 0. | 

Lemma 22. Suppose that we have a smooth vector field q = X(q) and a smooth map f in 
a neighborhood U C W 2 , satisfying condition If p = h(q), and h is a diffeomorphism 

between IA and h(U), then the induced vector field X(p) := (Dh) q X(q) and the map f = 
h o / o h^ 1 also satisfies condition 



det 



dVi dVi 

dx dy 

dV 2 dV 2 

dx dy 
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Proof. Indeed X(p) := (Dh) q X{q) = {Dh) h - Hp) X {h' 1 ^)). Hence 

X (/») = (Dh) h -i W(h -i m) X (h-Hhifih-Hp))))) = 
= (Dh) nh - Hq)) (Df) h . 1(p) X (h-Hp)) = 

= (Dh) f{h - Hp)) (Df) h - 1(p) (Dh-%] • [(Dh) h - Hp) X {h-\p))] 
= Df(p)-X(p). 



Lemma 23. Consider the planar vector field 
X{u,v) :■■ 



t 2 , 2\ ® , / 2 i 2\ ^ 

+u )—+ug(u +v ) — , 



(17) 



with g(0) ^ and / a differentiate map in a neighborhood of the origin, such that it leaves 
invariant the circles j r := {u 2 + v 2 = r 2 } and condition |7^p is satisfied, i.e. X(f(u,v)) = 
(D f) r UfV \X (u, v) . Then f is conjugated on each j r to a rotation with rotation number p{r) 
and limp(r) = 9, where cos9 ± isin(9 are the eigenvalues of Df\^ 0Q y 

Proof. By Theorem fl~8j and taking into account that the sets j r are also invariant under 



X, we obtain that on each j r , f is conjugated to a rotation and there exists t(\/u 2 + v 2 ) 
such that 

f(u, v) = ip (r(y/u 2 + v 2 ), (u, vj^j = ip (r(r), (u, v)) , 
where if is the flow of X. By taking polar coordinates it is not difficult to obtain that 



<p(t, (u,v)) 



cos 



(g(r 2 )t) - sin (g(r 2 )t) 
sin (g(r 2 )i) cos (g(r 2 )t) 



u 

V 



Hence on each set j r , f is indeed the rotation of angle a(r) = g{r 2 )r{r), which has 
rotation number p(r) = a(r)/(2ir). By using this fact and the differentiability of / at the 
origin we have 



f(u,v) 



cos a(r) — sina(r) 
sin a(r) cos a(r) 



u 

v 



cosa(O) — sina(O) 
sina(O) cosa(O) 



) 









+ 2 (u,v). 



Hence a(0) = 9 and therefore lim r ^o p( r ) = lim r ^o a(r)/(2ir) = 9/(2tt), as we wanted to 
prove. | 

Proof of Proposition[T9\ (i) From the above lemmas we know that the functions r(k, h, a) := 
r(p(k, h),a), and T(k, h, a) := T(p(k, h),a) are analytic functions in £; oc . Since by Theorem 
[18] we know that p(k,h,a) = r(k,h,a)/T(k,h,a), then the rotation number is analytic as 
well. 
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(ii) Consider the map F restricted on Q Since 
G(x, y, z) = x(x + l)z 2 + (x(x + 1) - y(y + l))z - y 3 - (1 + a + x)y 2 - (a + x)y, 
then equation G = is equivalent to z = z±(x,y) where 



y(y + l)-x(x + l)± y/A(x,y) 

z±{x,y) = —, — — r 

2x(x + 1) 

and A(x,y) = (y(y + 1) - x(x + l)) 2 + 4x(x + l)(y 3 + (1 + a + x)y 2 + (a + x)y). If x > and 
y > then > and Z-(x,y) < 0. Consequently the surface <5 can be described as: 

G = {(x,y,z+(x,y)) : x > 0, y > 0}. 

Hence, in a neighborhood of the fix point, F\g can be though as the planar map F(x,y) = 
(y, y)) in U = {(x, y) £ M 2 : x > 0, y > 0}. 

Clearly the map F(x,y) has (x c ,x c ) as a fix point and the matrix DF\, ^ , has the 
eigenvalues given by A = cos ± i sin where 



, l + x c \ /a-l + Vl + a 

^ = arccos — = arccos 



2x c J \ 2a 

Let X{x, y, z) be the vector field given by (|16p . Then the map -F(x, y) has an associated 
vector field 

d d 
X\ g =: X[x,y) = Xi(x,y,z+(x,y))— + X 2 (x, y, z+{x, y)) — , 

which is the restriction of the vector field (|16|) on ^. It can be checked that the vector field 
X has the point (x c , x c ) as a singular point, which is a non-degenerated center and satisfies 
X{F) = {DF)X. Via an analytic change of variables, X is conjugated with a vector field 
given in the normal form (|17p . say X^j. Through this conjugation we also obtain that F 
is conjugated with a new map Fjq (which by Lemma [22] satisfies condition (|14p with Xjv)- 
The maps JP and Fjv share the same eigenvalues at their respective fixed points. Using 
Lemma[23]we have that \\va.k^k c pp{k) = 9/(2tt), as we wanted to prove. 

(iii) The study of this case is similar to one of (ii), where here F and the surface Q are 
replaced by F 2 and Mh, respectively. 

A tedious computation shows that the characteristic polynomial of (DX) Ph is given by 
P(X) = A (A 2 + p(xh, a)) , where 

>2 



p(xh,a) 



(x h + 1) (2x h + a - 1) (-a + 3x h a - x h + 1 + 2x h 2 ) (x h 2 + x h + a - l) 



(a + x h ) 2 x h 2 (x h - l) 2 

Since x^ > 1 and a > 0, we have that p(xf l ;a) > 0, hence the eigenvalues of (DX) Ph 
are (which corresponds to the tangential direction of £), and a couple of conjugated 
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pure imaginary ones. By the implicit function theorem, in a neighbourhood of ph the set 
My' := Mh n {G > 0} is a differentiable manifold of dimension 2, invariant by X. Thus X 
restricted to induces a two dimensional vector field having a non-degenerated center 
at ph- At this point the proof follows in the same manner than in (ii). The computation of 
6h is straightforward. g 

The following result will be useful to study the odd periods of F. Although it seems 
natural that it is true for any a / 1, we have not been able to provide a general proof. 

Proposition 24. Consider a = a* := ^c-i)' 2 — 8-29590, where c = cos(27r/7). Then, there 
exists e > such that for any value of a satisfying \a — a*\ < e, the rotation number of F 
over the invariant curves {V\ = k} which foliate Q is not constant. 

Proof. First we prove the result for a = a* . We proceed by contradiction. If the set of 
rotation numbers were degenerated to a point, by Proposition [19] this value should be the 
value of the limit when we tend to the fix point, which is 



arccos g£ j _ 1 

2^ ~ 7' 

Indeed we have chosen a* to obtain this value. It gives the smallest denominator of all the 
rational numbers given by the expression (arccos ^ a ~ 1 "2^ 1 ~^ ) /(^ 7r )i where a € (0, oo). By 
Theorem [2] we also would have that F 7 restricted to Q would be the identity. 



On the other hand, take another point in Q, for instance q = (1, 1, x§<E£}M£+E^ ~ 
(1,1,3.20872). To prove that F 7 (q) ^ q it is convenient for the moment to consider F 
with a = ^d-i)' 1 anc ^ r = (1> 1> ^ 8d 2d~i d ~^ being d an unknown parameter. The equation 
that forces that the first components of F 7 (r) and r coincide is 

(51-272ci+540d 2 -464d 3 +144d 4 ) V8d 2 -12d+5- 1 W+72M- 1900d 2 +2504d 3 - 1664d 4 +448d 5 _ 
(38-168d+292(i 2 -240d 3 +80d 4 )v / 8(i 2 -12d+5-75+470d-1156cZ 2 +1400d 3 -832d 4 +192cZ 5 _ ' 

Working with the above equation we obtain that its solutions are included in the solutions 
of 

64(d - lf{d 2 - d/2 - 1/8) (d 2 - 3d/2 + 5/8)(128d 4 - 64<i 3 - 128d 2 + Wid - 19) = 0. 

Since the value d = c is not a solution, we have got that F 7 (q) ^ q, which is in contradiction 
with our initial assumption. Thus for a = a* we have proved that the set of rotation numbers 
on Q is not degenerated to a point. Recall that in Proposition [19] (i) it is proved that the 
rotation number varies continuously with respect to initial conditions and the parameter a. 
From this result we obtain that the set of all the rotation numbers over Q is not degenerated 
to a point for the values of a in some neighbourhood of a* , as we wanted to prove. Notice 
that when a = 1 this rotation number over Q, and also over + \ C, is reduced to the value 
1/8. | 
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5 On the set of periods of F. Proof of Theorem [3] 



In this section we prove Theorem [3] and its consequences (Corollary H] and Proposition [5]). 
Firstly, we present a constructive way for obtaining the denominators of the irreducible 
fractions which belong to a given interval. 

Theorem 25. Fix a real open interval I = (a,b) and denote by p\ = 2,p2 = 3,ps, . . . ,p n , . . . 
the set of all the prime numbers, ordered following the usual order. Associated to I we 
consider the following natural numbers: 

(i) The smallest prime number p m +i satisfying that p m +\ > max(3/(6 — a), 2), 

(it) Given any prime number p n , 1 < n < m, the smallest natural number s n such that 
Pn > 4/(6 -a). 

By using the above numbers, define the following finite subset of N : 

F sl>S2 ,..., Sm := {n £ N : n = p*fp% ■■■p^ with < U < - 1, % = 1, 2, . . . , m}. 
Then for any r G N \ i ? SliS2 ,.--.sm there exists and irreducible fraction q/r such that q/r £ /. 
Next result easily follows from the above Theorem: 

Corollary 26. Fix an open real interval (a,b). Following the notations of the above the- 
orem consider the number p := p^^p^ 2-1 • ■ 'P S m~^ ■ Then, for any r > p there exists an 
irreducible fraction q/r such that q/r £ (a,b). 



Proof of Theorem \25\ We prove the following two assertions: 

(a) If p is a prime number and p > p m +i then for any natural number k > 1 there exists 

an irreducible fraction of the form - — E I. 

kp 

(b) If pj is any prime number pi < p m +i and is the integer number given in the statement 
of the theorem, then for any natural number k > 1 there exists an irreducible fraction 

q 

kp\ 

Clearly the theorem follows from them. 

Let us prove the first one. From the fact that p m +i > 3/(6 — a), we have that if p is a 
prime number and p > p m +i then there exists an t such that 

a < < - < < b, (18) 

p p p 



of the form - € /. 
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where the three fractions are irreducible. Hence we have proved our assertion (a) for k = 1. 
Take now any k > 1. From the above inequalities we have that 

k£ — k ki-1 ki kl+l ki + k 
kp kp kp kp kp 

k£ — 1 kl+l 

Note that either — or — have to be irreducible because the factors of k never 

kp kp 

divides their numerators and if both were reducible the number p should divide both num- 
bers ki ± 1. Taking their difference we would have that p divides 2, a contradiction. Thus 
assertion (a) is proved. 

Let us prove assertion (b). Fix any prime number p = p n , smaller that p m +i and consider 
its associated number s = s n . From the inequality p^ n > 4/(6 — a) we have that 

a < < — < < < b, 

pS pS pS pS 

Note that either j + 1 or j have to be coprime with p hence taking £ either j or j + 1 
we have that 

£-1 £ £ + 1 , 
a < < — < < b, 

pS pS pS 

being the fraction £/p s irreducible, like in (|18p . When p > 2 we can argue as in the previous 

. , kl - 1 kt + 1 , 
case and assert that either — or — have to be irreducible, proving our result. 

kp s kp s 

When p = 2 we consider 

a < — < < < b. 

2 s 2 s 2 s 

Taking k > 2 we have that 

kj kj + 1 kj + 2 kj + k 

a< k¥ < ~klF~ < ~klF~ < ~1¥~ < ' 

kj ' + 1 kj-\-2 

and again one of the fractions — , — ; , has to be irreducible, as we wanted to prove. 

& k2 s ' k2 s ' ' F 



In the sequel we prove Theorem [31 Corollary U] and Proposition [5j 
Proof of Theorem^ For each a > 0, a ^ 1 and each x > 1 consider the function 



1 ({a-l){l-x) 
r(x) = — arccos 



2tt V 2 (ax + x 2 ) 

Recall that from Proposition [19] (hi), the function r{x) gives the limit of pp^ik) when k 
tend to Vi(p x ), where p x is the point on C given by (x, ,x). Observe that r(x) has a 
unique critical point which is a maximum (resp. a minimum) at x = x c = 1 + -^/l + a when 
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a > 1 (resp. < a < 1). Furthermore r(l) = 1/4 and lim a! _ > . 00 r(x) = 1/4. Now consider 
the value r(x c ) and denote it by p a : 

l ( (l-a)vTT^ \ , 1Q , 

Pa = — arccos — . , . (19) 

Take a > 1 and a number p* G (1/4, p a ) (the case a < 1 and p* G (p a , 1/4) can be studied 
in a similar way). Let us see that there is a continuum of initial conditions in {G > 0} \ C 
such that their rotation number is p* (and notice that by using expression ([6]), the images 
by F of these initial conditions satisfy the same property and are in {G < 0} \C). For e > 
small enough, there are two periodic points of F in £, say p^ = (x^, ,x^), such that 
r (p ± ) = P*±£- By Proposition[l9](i) there exist initial conditions r^ G {Vi = Vi(p =1: )}n{G > 
0} such that their respective rotation numbers, g^ satisfy p* — 2e < g~ < p* < g + < 
p* + 2e. Joining r~ and r + by a continuous path T C {G > 0} \ C, and by using again the 
continuous dependence of the rotation number with respect the initial conditions, we obtain 
the existence of a point r G V such that its rotation number is exactly p* . By Theorem [2] 
the same happens with all the points in {G > 0} of {Vi = V\(r)} n {V2 = ^(r)} = S 1 , as 
we wanted to prove. g 

Proof of Corollary [7J (i) By using Theorem [3] and Corollary [26] the result follows, 
(ii) Observe that the function p a given in (fl~9|) is an increasing function such that 

v vr - 2 arcsin(l/8) 1 
hm p a = ~ 0.23005 and lrm p a = -. 

a^0+ 47T 3 

Therefore, by using again Theorem [3j for each number in ((tt — 2 arcsin (l/8))/(47r), 1/4) 
there exists some a G (0, 1) and some initial condition outside Q with this rotation number 
for F 2 . Similarly, for each number in [1/4,1/3) there exist some a > 1 and some initial 
condition, also outside Q, with this rotation number. In particular, for all the irreducible 
rational numbers p/q with the property 

-/r — 2 arcsin (1/8) p 1 
4vr < q < 3 

we can find a value of a such that F 2 has continua of periodic orbits of period q. 

(Hi) Setting a = (tt — 2 arcsin (1/8) )/(4ir), b = 1/3 and using the notation introduced 
in Theorem 125} we have that, m = 10, pn = 31 and p\ = 2 (with s\ = 5), P2 = 3 (with 
S2 = 4), p 3 = 5 (with s 3 = 3), p 4 = 7, p 5 = 11, p 6 = 13, p 7 = 17, p 8 = 19, p 9 = 23 and 
P10 = 29 (where Si = 2 for i G {4, . . . , 10}). From Theorem 125} we have that for all q G N, 
such that q > q := 2 4 • 3 3 • 5 2 • 7 • 11 • 13 • 17 • 19 • 23 • 29 = 2 329 089 562 800 there exists 
some a > and some g-periodic orbit for F 2 . It is now easy to develop a finite algorithm 
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in order to find which irreducible fractions p/q with q < qo are in I rot . Implementing this 
algorithm we get that there appear irreducible fractions with all the denominators except 
1,2,3,5,6,8,9,12,14 and 20. Doubling these numbers, and taking into account that C is 



Proof of Proposition [5l (i) This result is a direct consequence of expression (|6|) . 

(ii) Fix a value of o of the ones given in Proposition [MJ By this value we know that 
the set of all rotation numbers of all the points of Q contains an open interval. By applying 
Corollary 1261 to this interval the result follows. 

(iii) Similarly that in the proof of (i) of Theorem [3l for each a > 0, we consider the 
function s(a) = ^ arccos ^ a ~ 1 ^ 1 — ^ . Recall that from Proposition [19] (iii), this function 
gives the limit of the rotation numbers over Q when we approach to the fix point. The 
range of this function when a > is J ro t • Taking into account the continuity of the rotation 
number with respect to initial conditions and the parameter a, and arguing as in the last 
part of proof of Theorem [3j the result follows. g 

6 Some numerical results 

In this section we present some numerical explorations which lead us to establish the open 
questions stated in Section fOl 

The following tables of rotation numbers have been obtained using the relation (|15p . 
by numerical integration of the vector field (|16|) using a 7-8^ order Runge-Kutta method. 
Table 1 has been obtained taking a = 3, and gives the rotation number associated to 
the orbit passing through some points of the surface Q. These points have been taken by 
considering the following path over Q: 



where z(x(t),y(t)) is one of the two branches of solutions of equation G(x(t),y(t), z) = 0. 
Recall that r(p) := Tp(p) is given by the relation ip(Tp(p),p) = F(p), T(p) is the period 
of the periodic orbit of (|16p passing through p and the rotation number of F at the orbit 
starting at p is pf(p) = t f(p)/T(j>). Similarly we can define r F 2(p) and pp2(p). Note also 
that when both numbers have sense pp2 {p) = 2pf(p)- In general we observe that the function 
t — > p(x(t),y(t), z(t)) seems to be decreasing. 



full of two periodic points of F, we obtain (iii). 



I 




(20) 



30 



t 


Point p 


T \P) 


tf(P) 


Pf(P) 


P F 2 (p) 





(3,3,3) 






0.13386 


0.26772 


1 


(4,7/3,1.62395) 


0.41781 


0.05586 


0.13369 


0.26737 


2 


(5,2, 1.06969) 


0.36063 


0.04810 


0.13337 


0.26674 


3 


(6,9/5,0.78049) 


0.30622 


0.04074 


0.13305 


0.26610 


4 


(7,5/3,0.60637) 


0.26009 


0.03452 


0.13274 


0.26549 


5 


(8,11/7,0.49153) 


0.22226 


0.02944 


0.13247 


0.26494 


6 


(9,3/2,0.41083) 


0.19148 


0.02532 


0.13223 


0.26446 


7 


(10, 13/9,0.35143) 


0.16635 


0.02196 


0.13201 


0.26402 


8 


(11,7/5,0.30610) 


0.14570 


0.01921 


0.13182 


0.26364 


9 


(12,15/11,0.27051) 


0.12860 


0.01693 


0.13164 


0.26328 


10 


(13,4/3,0.24191) 


0.11430 


0.01503 


0.13149 


0.26298 


11 


(14,17/13,0.21849) 


0.10225 


0.01343 


0.13134 


0.26269 


12 


(15,9/7,0.19899) 


0.09202 


0.01207 


0.13122 


0.26244 


13 


(16,19/15,0.18253) 


0.08325 


0.01091 


0.13110 


0.26220 



Table 1. Rotation number on Q for a = 3. 



Note that the results of Table 1 also give light to know which are the odd periods for F for 
a given value of a. For instance when o = 3 it seems clear that 0.1333 . . . = 2/15 is one of the 
rotation numbers reached by F over Q. Hence for this value of a, F must have periodic points 
of period 15. By applying a three dimensional Newton method to the system F 15 (x, y, z) = 
(x,y,z) we have obtained the approximated solution r ~ (2.00557,5.20647,9.89389). Note 
that \F 15 (r) — r\i < 3.8 x 10~ 5 , where as usual \(x, y, z)\\ = \x\ + \y\ + \z\. Indeed there should 
exist infinitely many 15-periodic points, given by all the orbits starting at the periodic orbit 
of (|16p passing through a given 15-periodic point. We also have checked that 

1^15x10"^ _ < 2 8 x 1Q n-6 for n = 1)2 ,3,4. 

If r were a true 15-periodic point the above values should have been zero, but as we have 
already noticed in Remark [6] the dynamical system generated by F has sensible dependence 
with respect to initial conditions. 

Table 2 is again obtained taking a = 3. Now the rotation number of F 2 is computed for 
some points in the curve given by 

p(t) = {x(t),y(t), z{t)) = (ix c + 9) • (1 - t) +tx 1: ^0l , z {x(t),y(t)) \ , t G [0, 1], 

(21) 

where z(t) is a fixed branch of the two branches of solutions of equation (|7|) and x(l) = x\ 
is the first coordinate of the two periodic point of F, p(l) ~ (1.11929,34.53097,1.11929). 
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This curve joints the point p(0) = (12, 15/11, 0.27050 . . .) E Q with p(l) running over the 
level surface 

= 28561 (19 + 3V89)(197 + VS) 

43560 -79 + 13 v/89 





Point p 


1 \P) 


T F 2 (p) 


Pf 2 (p) 





(12, 1.36364,0.27051) 


0.12860 


0.03386 


0.26328 


0.1 


/iaaiiao 1 /lAotrir a o /i to t\ 

(10.91193, 1.40355, 0.24737) 


A 1 O O C T 

0.12857 


A AO O O tr 

0.03385 


0.26327 


0.2 


(9.82386,1.45332,0.22528) 


0.12848 


0.03382 


0.26325 


0.3 


(8.73579,1.51708,0.20426) 


0.12829 


0.03376 


0.26320 


0.4 


(7.64772,1.60171,0.18437) 


0.12796 


0.03367 


0.26310 


0.5 


(6.55965,1.71947,0.16576) 


0.12743 


0.03351 


0.26295 


0.6 


(5.47158,1.89454,0.14879) 


0.12655 


0.03324 


0.26270 


0.7 


(4.38350,2.18221,0.13432) 


0.12501 


0.03279 


0.26226 


0.8 


(3.29543,2.74259,0.12498) 


0.12205 


0.03191 


0.26145 


0.9 


(2.20736,4.31300,0.13351) 


0.11498 


0.02985 


0.25962 


0.95 


(1.66333,7.03020,0.17364) 


0.10648 


0.02744 


0.25768 


0.99 


(1.22810,18.53618,0.43340) 


0.09080 


0.02314 


0.25484 


0.999 


(1.13017,31.72824,0.95253) 


0.08589 


0.02183 


0.25414 


0.9999 


(1.12038,34.22789,1.09981) 


0.08575 


0.02179 


0.25412 


1 


(1.11929,34.53097,1.11929) 






0.25412 



Table 2. Rotation number on {Vi = k* ~ 146.70452} when a = 3. 



Table 3 has been obtained by repeating the first experiment but taking a = 7/9. So it 
gives the rotation number associated to the orbit passing through some points of the surface 
Q, by considering the path given by (|20p . Notice that in this case the rotation number seems 
an increasing function of t. 
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t 


Point p 


T KP) 


tf{p) 


Pf(P) 


P F 2 (p) 





(7/3,7/3,7/3) 






0.12338 


0.24676 


1 


(10/3,37/21,1.11361) 


0.48969 


0.06043 


0.12340 


0.24681 


2 


(13/3,23/15,0.71973) 


0.39978 


0.04935 


0.12345 


0.24690 


3 


(16/3,55/39,0.52965) 


0.32588 


0.04025 


0.12350 


0.24700 


4 


(19/3,4/3,0.41853) 


0.26908 


0.03324 


0.12354 


0.24708 


5 


(22/3,73/57,0.34583) 


0.22552 


0.02787 


0.12358 


0.24716 


6 


(25/3,41/33,0.29462) 


0.19171 


0.02370 


0.12361 


0.24722 


7 


(28/3,91/75,0.25662) 


0.16502 


0.02040 


0.12364 


0.24729 


8 


(31/3,25/21,0.22731) 


0.14363 


0.01776 


0.12367 


0.24734 


9 


(34/3,109/93,0.20401) 


0.12622 


0.01561 


0.12369 


0.24739 


10 


(37/3,59/51,0.18506) 


0.11187 


0.01384 


0.12372 


0.24744 


11 


(40/3,127/111,0.16933) 


0.09990 


0.01236 


0.12374 


0.24748 


12 


(43/3,17/15,0.15607) 


0.08980 


0.01111 


0.12376 


0.24752 



Table 3. Rotation number on Q for a = 7/9. 
Finally, Table 4 is obtained taking a = 7/9, considering the path (|2ip . which runs 
from the point p(0) = (34/3, 109/93, -296730/604469 + 1/1813407^1587984839746) ~ 
(11.33333, 1.17204, 0.20401) e Q to the fixed point of F 2 , p(l) ~ (1.04794,38.08255,1.04794), 
over the level surface 



101 (923217+V1587984839746) (69592724+3 V1587984839746) 
■ — 272452149 -890190+ V1587984839746 



t 


Point p 


Tip) 


T F 2 (p) 


p F 2 (p) 





(11.33333,1.17204,0.20402) 


0.12622 


0.03123 


0.24738 


0.1 


(10.30479,1.19106,0.18573) 


0.12620 


0.03122 


0.24739 


0.2 


(9.27625,1.21480,0.16809) 


0.12612 


0.03120 


0.24740 


0.3 


(8.24772,1.24529,0.15108) 


0.12596 


0.03116 


0.24742 


0.4 


(7.21918,1.28585,0.13473) 


0.12569 


0.03110 


0.24744 


0.5 


(6.19064,1.34250,0.11911) 


0.12526 


0.03010 


0.24748 


0.6 


(5.16210,1.42714,0.10435) 


0.12455 


0.03083 


0.24755 


0.7 


(4.13356,1.56733,0.09081) 


0.12335 


0.03055 


0.24765 


0.8 


(3.10502,1.84454,0.07963) 


0.12107 


0.03001 


0.24784 


0.9 


(2.07648,2.65147,0.07637) 


0.11551 


0.02867 


0.24824 


0.95 


(1.56221,4.16212,0.09117) 


0.10807 


0.02687 


0.24866 


0.99 


(1.15080,12.78937,0.23101) 


0.08942 


0.02229 


0.24932 


0.999 


(1.05822,31.53212,0.74344) 


0.07887 


0.01968 


0.24955 


0.9999 


(1.04897,37.30369,1.00598) 


0.07829 


0.01954 


0.24956 


1 


(1.04794,38.08255,1.04794) 






0.24957 



Table 4. Rotation number on {V\ = k ~ 0.24956} when a = 7/9. 
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Appendices 



A Proof of Lemma [12 

To describe the foliation of Q + , induced by A(x,y;a,k) = obtained for a fixed value of 
a > 0, and varying k > k c , we solve the quadratic equation (with respect k): A(x, y; a, k) = 
x 2 y 2 k 2 + pi(x, y;a)k + po(x, y\ a) = 0. Thus the curve A(x, y; a, k) = in Q + can also be 
described by two equations 

, (x + y + a + 1 ± 2 + y + a) (x + 1) (y + 1) 

A; = m±(x, y; a) = - . 

xy 

We make the following claims: 

Claim 1: For any fixed k >k c the following statements hold 

(i) If a > 1, there exist two values xn~ < xi,k such that equation (with unknown y) 

k = m-(x,y;a), (22) 

has two solutions if x € {xiki x 2k)j one solution if x = Xik i = 1,2, and none solution if 
% ^ [ x \ki x 2k\- This means that varying x > 0, equation \22\) describes an oval 

(ii) If a < 1, then there exist a value x^ > 1 — a such that equation \22\) has two solutions 
if x < Xk, one solution if x = x^, and none solution if x > x/.. This means that varying 
x > 0, equation h22\) describes a curve consisting of a point (xk,yk) an d (from right to 
left) two positive branches yi(x) < y2{x) defined only for x G (0,Xfc). A more accurate 
analysis will show that these two branches meet at the point (x,y) = (0, 1 — a). Therefore 
they describe an oval, namely Cfc- 

Claim 2: (i) For k > k c , there exist two values x\^ < x c < X2,fc such that equation 

k = m + (x,y;a), (23) 

has two solutions if x £ (xi t k, x 2,k), one solution if x = x^k i= 1,2, and none solution if 
% ^ [ x i ki x 2k\- This means that varying x > 0, equation I123\) describes an oval 7^. 

(ii) The equation k c = m + (x,y; a) has a unique solution if x = x c and none solution if 

Since m + (x,y; a) > m_(x,y;a) it is easy to see that each oval 7& surrounds the corre- 
sponding oval Cfc- From this fact and the above claims the proof of the lemma follows. 
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Before proving both claims we establish some common facts. We fixi>0 and we use 
the following notation: 



dm+ ,_ , , 1 . 
— — (x,y;a) = -{x + 1) 
dy 



±f(x,y;a) + g(x,y;a) v / h(x,y;a) 
xy 2 ^h(x,y;a) 



where f(x, y; a) = —y 2 + y + 2x + 2a, g(x, y; a) = — y 2 + x + a + 1 and /i(x, y; a) = x + y + a. 

The solutions in Q + of — — (x,y;a) = are described by 

dy 

(f 2 ~ 9 2 h) {x,y;a) = {y + x + a- l)(-y 2 -y+x+a- l)(-y 2 + y + x + a). (24) 

So this equation gives the local extrema of y —* m±(S, y; a). 
It can be easily proved that 

lim m±(x, y; a) = +oo and lim m±(x, y; a) = +oo, (25) 

for all x > 0. 

Let us now proceed with the proof of both claims. 

Proof of Claim 1. (i) If a > 1, and taking into account that x > 0, it is easy to see that 
— — (x, y; a) = if and only if 

Q(x,y) = -y 2 -y + x + a- 1 = (26) 

which has the unique positive solution y = y m i n (x) = (-1 + V-3 + 4a + 4x)/2. 

We point out that this solution is well defined for all x > and a > 1, and that 
ymin(x) > 0. Hence, taking into account equation (|25|). the function y — > m_(x, y; a) takes a 
minimum at the point y m in(x). So we have that for each x > the functions y — > m_ (x, y; a) 
are decreasing in the interval y £E (0, ymin 

) and increasing in y £ (y m in> +oo). 
Now we study the function x — > m_(x,y m i n (x);a). We have the following facts: 

(I) lim m_(x,y min (x);a) = +oo, 

(II) It is easy to see that at infinity m_(x, y m in(x); a) ~ x, thus 

lim m_(x,y min (x);a) = +oo. 



x— >+oo 



Qffi 

(III) Since the only positive solutions of — - — (x, y; a) = are given by the solutions of the 

ox 

equation Q(y, x) = —x 2 — x + y + a— 1 = 0, where Q(x, y) is defined in ([26]) . and Q(x, y) = 

Q(y,x) = if and only if x = y = x* := y/a — 1, we obtain that - ~ (x, y m i n (x); a) / 

ox 

Qfn 

for x ^ x*, and — — (x*,x*:a) = 0. This means that x* is the unique critical point of the 

ox 

function x — ► m_(x, y m j n (x); a) when x > 0. 
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(IV) It can be easily checked that fc* = m_(x*, x*; a) = (1 + V a — I) 2 an d ^* < Re- 
collecting the results in items (I)-(IV) we obtain that the function x — > m_ (x, y m in(x)',a), 
which gives the minimum values of each function y — > m_(x, y; a), has a unique minimum 
at x = x*, decreases from +00 to fc* for x £ (0,2;*), and increases from to +00 for 
x G (x*, +00). 

This proves that for any fixed k > k c there always exist two solutions X{k , i = 1,2 of 
equation A: = m_(x, y m i n (x); a), see also Figure 3. 

For these two values, the minimum of the functions y — > m_ (xj y; a) is fe. This implies 
that equation (|22|) only has one solution for x = xn- , i = 1,2. Now observe that for all 
x £ (xi^, X2,fc)! since the minimum values of the functions y — > m_(x, y; a) are below k, we 
can conclude that for these values of x equation (|22p only has two solutions. Finally for 
x ^ 0^1 fc> ^2,^)) since the minimum values of the functions y — > m(x,y;a) are greater than 
k, equation (|22p has no solutions. In summary, equation (|22|) describes one and only one 
oval Cfc- 

a>\ a < 1 




2/mmM ^ «imm(l) = 1 - a - X >" 

a > 1 or a < 1 with x > 1 — a a < 1 with x < 1 — a 

Figure 3. Plot of some functions involved in the proof of Lemma [T2l 
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(ii) If a < 1, — — (x, y; a) = at the curve y = y^nix), described by the unique positive 

oy 

root of (y + x + a — 1)(— y 2 — y + x + a — 1). Thus 



1 — a — x if x < 1 — a, 

(-1 + V-3 + 4a + 4x)/2 if x > 1 - a, 



which is well defined for all x > 0. Taking into account equation (|25h . y m i n (x) gives a 
minimum for y — ► m_(x, y; a). 

Now, we look at the function x — ► m_(x, y m i n (x); a). First notice that for x < 1 — a, 
m_(x, y m in(^); «) = m-(x, x + a — 1; a) = 0. In the region x > 1 — a, we have 



"i_(x,y min (x);a) 



n(x)(x + !)(! + \/-3 + 4a + 4x) 2 



8x(x — 1 + a) 



where n(x) = f(s(x)) with s(x) = 2(x + a) + \/-3 + 4x + 4a, and [(s) = s + 1 + 2\/2V 's - 1. 

Observe that x > 1 — a if and only if s > 3, and that the function is monotonic 
increasing from to +oo for s > 3. Hence for x > 1 — a, we have m_(x, y m i n (x); a) is 
monotonic increasing from to +oo. 

Therefore, for x < Xk, the minimum of the functions y — > m_(x, y; a) is always below k, 
hence we can conclude that for these values of x equation (|22p has exactly two solutions, 
giving rise to two positive branches y\{x) and yi (x). Since A(0, y; a, k) = (y + l) 2 (a — 1 + y) 2 
we can conclude that these two branches meet at the point (0, 1 — a), which is an order two 
contact point of with {x = 0}. The symmetry of A with respect the line {y = x}, gives 
the other contact point (1 — a, 0). 

For x = Xfc the minimum of the functions y — > m_(x,y;a) is k. This implies that 
equation (|22p has only one solution, giving the point where the previous mentioned two 
branches meet. 

Finally, observe that for all x > Xk, since the minimum values of the functions y — * 
m_(x,y;a) are over k, equation (|22[) has no solutions. 

In summary equation (|22p describes one and only one oval which has two contact 
points with the boundary of Q + at (0, 1 — a) and (1 — a, 0). This ends the proof of Claim 1. 

Proof of Claim 2. We have to study the function m+. It is easy to see that y — > m + (x, y; a) 
has a unique minimum at y = y m i n (x) = 1 +Vi+4a+ 4g ^ ^ e unique positive solution of the 
equation C(x,y) := — y 2 + y + x + a = 0, and that m + (x, y; a) is decreasing from infinity to 
m + (x, ymin(x); a) for y < y min (x) and increasing to infinity for y > y min (x). 
We need now to study the function x — > m + (x, ymm(x); a). Observe that: 

(I) It holds that 

(II) Since lim y m \ n (x) = (1 + y/l + 4a)/2 > 0, we get lim m + (x, y m \ n (x); a) = +oo. 

a;->0+ x^0+ 
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(Ill) It is easy to check that at infinity m + (x, y m i n (x); a) ~ x, thus 

lim m + (x,y min (x);a) = +00. 



x— >+oo 



(IV) Since the only positive solutions of — — (x,y;a) = are given by the equation 

ox 

C(y,x) = —x +x + y + a = 0, and C(x,y) = C(y,x) =0 if and only if x = y = x c = 
1 + \/T+~a, we have that (x, y m in(^); a) / for x / x c , and (x c , y m in(a; c ); a) = 0. 

Collecting the information summarized in (I)-(IV), we obtain that the function x — * 
m + (x,y Ta i n (x);a), which gives the minimum values of each function y — > m + (x,y; a), has 
a unique minimum at x = x c , and decreases from +00 to k c for x £ (0, x c ), and in- 
creases from k c to +00 for x £ (x c ,+oo). A simple computation, omitted here, shows that 
m + (x c , y min (x c ); a) = m + (x c , x c ; a) = k c . 

This proves that for any fixed k > k c there exists only two solutions Xj&, i = 1,2 of 
equation 

k = m + (x,y min (x);a), (27) 

such that x c E (^i,fc,^2,fc)) see again Figure 3. For these two values, the minimum of the 
functions y — * m+(xi t k,y,a) is k. This means that equation (|23p has only one solution for 
x = i = 1,2. We note that for all x G (^ifej^2fc)) since the minimum values of the 
functions y — > m + (x, y; a) are below fc, we can conclude that for these values of x equation 
(|23|) has only two solutions. Finally for x ^ {x\ ^, Xi k), since the minimum values of the 
functions y — * m + (x,y; a) are greater than k, equation (|23p has no solutions. In summary 
equation (f23|) describes one and only one oval 7^. This ends the proof of (i). 

The above analysis of x — > m + (x, y m i n (x); a) shows that k c = m(x,y;a) if and only if 
x = y = x c . So the branch of Tk c described by equation ([23]) collapses to the point (x c , x c ). 
This ends the proof of (i), and so the proof of Claim 2. g 



B Proof of Lemma [13 

To describe the foliation of Q + , induced by A(x,y;a,h) = obtained for a fixed value of 
a > 0, and varying /i > /i c , we can rewrite A(x, y; a, h) = x 2 y 2 /i 2 +pi(x, y; a) h +po(x, y; a) 
where 

Po(x, y; a) = x 2 y 4 + (2x 3 + 2x 2 + (-2a + 2) x) y 3 

+ (x 4 + 2x 3 + (-4a + 5) x 2 + (-4a + 4) x + a 2 - 2a + l) y 2 

+ ((-2a + 2) x 3 + (-4a + 4) x 2 + (4 + 2a 2 - 6a) x + 2 + 2a 2 - 4a) y 

+ (a 2 - 2a + l) x 2 + (2 + 2a 2 - 4a) x + a 2 - 2a + 1, and 

pi(x, y, a) = ((-2x 2 - 4x))y 3 + (-2x 3 - 10x 2 + (-2a - 6)x)y 2 
+ (-4x 3 + (-2a - 6) x 2 + (-2 - 2a) x)y. 
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Thus the curve A(x, y; a,h) = in Q + can also be described by two functions 

(yx + 2x + 2y + a + 1 ± 2 v / d(x,y;a)) (1 + x + y) 

h = m±(x,y;a) = ± '- , 

yx 

where 

d(x, y; a) = x 2 y + xy 2 + x 2 + y 2 + (a + 2)xy + (a + l)x + (a + l)y + a. 
As in the previous appendix, to prove the lemma we make two claims: 

Claim 1: For h > h c and for all fixed x > 0, we will see that there exist two solutions 
yi(x) < y2(x) of the equation 

h = m_(x, y; a), (28) 

that give rise to the two branches of yi(x) and y2{x) given in the statement of the 
lemma. Moreover, lim ydx) = +oo and lim ydx) = + . 

Claim 2: (i) For h > h c , if we consider the equation 

h = m + (x,y;a), (29) 

then there exist two values x\^ < x^h such that equation \29\) has two solutions if x S 
(xi h, x% h)i one solution if x = Xih i = 1,2, and none solution if x ^ [xih, x 2h\- This 
means that varying x > 0, equation l2§\) describes the oval ^h- 

(ii) The equation h c = m + (x,y;a) has a unique solution if x = x c and none solution if 
x / x c . 

Note that, since m + (x, y; a) > m^(x, y; a), the solutions of the equation (|2"9"|) (whenever 
they exist) are contained in the interval (yi(x), y2(x)) defined by equation ([28]) . This implies 
that the oval 7/1 is contained between the two branches yi(x) and 2/2 (^)- So, by using the 
above two claims, the lemma follows. 

Before giving the proof of the above claims we establish some common facts. We fix 
x > and we use the following notation: 

dm± Tf(x, y;a) + g(x,y;a)^h(x, y; a) 

— — {x,y;a) = = , 

Oy xy 2 yh{x,y;a) 

where 

f(x, y; a) = {-2x -2)y 3 + (-x 2 + (-2 - a) x - 1 - a) y 2 
+ (x 3 + (3 + a) x 2 + (2a + 3) x + 1 + a) y 
+2x 3 + (2a + 4) x 2 + (2 + 4a) x + 2a, 

g(x, y\ a) = (x + 2) y 2 — 2x 2 + (—a — 3) x — 1 — a, 

h(x,y;a)= (y + 1) (x + 1) (x + a + y) . 
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So — — (x, y; a) = if and only if 
dy 

f 2 (x, y; a) - g 2 {x, y; a) h(x, y; a) = (x + 1) (-y 3 - y 2 + (1 + x)y + (1 + x){a + x)) ■ 
■{xy - a + 1) (xy 2 + (1 + x)(a - 1 + x)y + (1 + x)(a - 1)) = 0. 

This equation gives the local extrema of y — > m__(x, y; a) and y — > m+(x, y; a). 
It is not difficult to see that 

lim m±(x, y; a) = +oo and lim m±(x, y; a) = +oo, (30) 

for all x > 0. 

Proof of Claim 1. If a > 1, taking into account that x > 0, it is easy to see that 
Qrn 

— — (x,y;a) = if and only if yo = (a — l)/x. Furthermore m—(x,(a — l)/x;a) = 
dy 

for all x > 0, which is a minimum of y — > m__(x,y;a). Taking into account equation (|30p . 
we have that for each x > the functions y — > m_(x,y;a) are decreasing from +oo to 
for y G (0, yo) and increasing from to +oo if y € (yo, +oo). This proves that in this case 
equation (I28p always has two solutions for any h > 0, in particular for any h > /i c . 

If < a < 1, — — (x, y; a) = if and only if y = y q (x), where y q (x) is the only positive 
solution of the quadratic equation xy 2 + (1 + x)(a — 1 + x)y + (1 + x)(a — 1) = which 



taking into account (|30|) gives a minimum of y — > m_(x,y;a). To see that equation (|28p 
always has two solutions for h > /i c , we only have to see that m-(x,y q (x); a) < h c for all 
x > 0. Since m_(x, y; a) = 1 — a if and only if xy 2 + (1 + x)(a — 1 + x)y + (1 + x)(o — 1) = 0, 
we have that m-{x, y q (x); a) = 1 — a. On the other hand it is not difficult to check that 
1 — a < h c . So the first part of the claim is proved. 

To end the proof it remains to see that lim yAx) = +oo and lim yAx) = + , where 
yi(x) and ya(x) are the two branches of A = 0. We observe that each curve A(x, y; a, h) = 
is symmetric with respect the axis y = x. So it is equivalent to see that lim y«(x) = , 

x— >+oo 

i = 1, 2 or lim y^(x) = +oo i = 1, 2. So we will prove the first equality. To do this we will 
study if there arrive any branch of I\ to the "infinity line" in the projective space PR 2 in 
the direction y = 0. 

For each affine curve A(x,y;a,h) = we can consider the projectivized curve in PR 2 
(given in homogeneous coordinates (x,y,u)) Th = {A(x,y,u;a,h) = 0}, where 
A(x, y, u; a, h) = (a 2 - 2a + l) n 6 + (2 + 2a 2 - 4a) xu 5 + (a 2 - 2a + l) x 2 n 4 + 
(2 + 2a 2 - 4a) yu 5 + (4 + 2a 2 - 6a - 2h - 2ka) xyu 4 + (-4a - 2ka + 4 - 6h) x 2 yu 3 + 
(2 - 4h - 2a) x 3 yu 2 + (a 2 - 2a + l) y 2 n 4 + (-4a - 2ka + 4 - 6h) xy 2 u 3 
+ (h 2 - 4a - lOh + 5) x 2 y 2 u 2 + (-2h + 2) x 3 y 2 u + y 2 x 4 + (2-4h- 2a) xy 3 u 2 
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+ (— 2h + 2) x 2 y 3 u + 2y 3 x 3 + y 4 x 2 . In the local chart {x ^ 0} this curve is given by (just 
taking x = 1) = {A(l, y, u; a, h) = 0}, where: 

A(l, y, u; a, h) = (a 2 - 2a + l) u 6 + (2 + 2a 2 - 4a) u 5 + (a 2 - 2a + l) u 4 + 

(2 + 2a 2 - 4a) yu 5 + (4 + 2a 2 - 6a - 2h - 2ka) yu 4 + (-4a - 2ka + 4 - 6h) yu 3 + 

(2 -Ah- 2a) yu 2 + (a 2 - 2a + l) y 2 u 4 + (-4a - 2ka + 4 - 6h) y 2 u 3 + 

(h 2 - 4a - lOh + 5) y 2 u 2 + (-2h + 2) y 2 u+y 2 +(2 -Ah- 2a) y 3 u 2 +(-2h + 2) y 3 u+2y 3 +y 4 . 

We want to prove is that there are two "affine" branches of f\ arriving at the point of 
the infinity line with coordinates (y, u) = (0, 0) € r\. Since A(l, y, u; a, h) = y 2 + 2y 3 + (2 — 
2h)y 2 u + 2(1 — a — 2h)yu 2 + 0((x, y) 4 ), we need to perform the blow-up: (y, u) = (v u, u), 
which after removing the factor u 2 transforms Th into 
f* = { [ (a 2 - 2 a + 1) v 2 + (2 + 2 a 2 - A a) v + a 2 - 2 a + l] u 4 

+ [(2 - 4 h - 2 a) v 3 + (-4 a - 2 /ea + 4 - 6 a) w 2 + (4 + 2 a 2 - 6 a - 2 h - 2 ka) v + 2+ 
2a 2 -4a] n 3 + [v 4 + (-2 h + 2) v 3 + (h 2 - A a - 10 h + 5) u 2 + (-4 a - 2 fea + 4 - 6 a) v+ 
a 2 - 2 a + l] u 2 + [2 v 3 + (-2 a + 2) v 2 + (2 - 4 a - 2 a) v] u + v 2 = 0}. 

The intersection with {v = 0} are the points (v,u) = (0,0), and (v,u) = (0,-1). 
But the last point is not interesting for us since, the affine region {(x,y) : x > 0,y > 0} in 
this local coordinates corresponds with {(v,u) : v > 0,u > 0}. The directions of approach 
of to (v, u) = (0, 0) are given by: u = X± v, where 

(2h + a - 1) ± 2 v / /i(/i + a - 1) 

It is easy to check that for a ^ 1, both \± are positive, therefore there exist two branches of 
T* h arriving at the singular point (0,0) in {(v,u) : v > 0, u > 0}. As only the two branches 
of I\ described by y«(x), i = 1, 2 are defined when x — > +oo these ones are the two branches 
described by the blow-up procedure. This ends the proof of the claim. 

Proof of Claim 2. It is easy to see that y — ► m+(S, y; a) has a unique minimum at y = y c (x), 
where, by Descartes' rule, y c (x) is the only positive solution of the equation C(x,y) := 
—y 3 — y 2 + (1 + x)y + (1 + x){a + x) = 0, and its decreasing at y < y c (x) at increasing for 

y > y c {x~). 

Now we state three facts concerning the curve y = y c (x), which are relevant to the study 
of x — > m + (x, y c (x); a). 

(a) An straightforward computation shows that y c (x c ) = x c . 

(b) Once again, applying Descartes's Rule on the cubic C(x, y) = 0, we have that lim y c (x) > 
0, which implies that lim m + (x, y c (x); a) = +oo. 
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(c) A detailed analysis of the asymptotic expansion of y c (x) at infinity gives that y c (x) ~ 
(•^216/6) V?, hence lim y c (x) = +00, and as a consequence lim m+(x, y c (x); a) = 

X—r + OO X— > + 00 

+ OO. 

Since the only positive solutions of — — (x, y; a) are given by the cubic equation C(y, x) = 

—x — x 2 + (1 + y)x + (1 + y)(a + y) = 0, and C(x, y) = C(y, x) = if and only if x = x c 

and y = x c , we get that — — (x,y c (x);a) ^ for x ^= x c , and — — (x x ,y c (x c );a) = 0. 

ox ox 
This means that the function x — ► m + (x,y c (x); a), which gives the minimum values of 

each function y — > m+(x, y; a), has a unique minimum at x = x c , decreases from +00 to 

h c for x G (0, x c ), and increases from /i c to +00 for x G (x c ,+oo). A simple computation 

shows that m + (x c , y c (x c ); a) = m + (x c ,x c ;a) = /i c . 

The above results prove (see Figure 5), that for any fixed h > h c there exists only two 

solutions Xih, £ = 1, 2 of equation 



h = m + (x,y c (x);a), 



(31) 



such that x c € {xih,%2h)- For these two values, the minimum of the functions y — * 
m{xih,y; a) is h. This means that equation (|29p has only one solution for x = Xi t h i = 1, 2. 
Now we observe that for all x E (xi h, ^2 ft)> since the minimum values of the functions 
y — > m(x, y; a) are below /i, we can conclude that for these values of x equation (|29p only 
has two solutions. Finally for x ^ (xihi x 2h)i since the minimum values of the functions 
y — > m(x, y; a) are greater than h, equation (|29p has no solutions. In summary equation 
(|29p describes one and only one oval 7^. 





Figure 5. Plot of some functions involved in the proof of Lemma [T3l 
The above analysis of x — > m + (x, y c (x); a) shows that /i c = m(x,y;a) if and only if 
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x = y = x c . So the branch of T/ lc described by equation (f29|) is only the point (x c , x c ). This 
ends the proof of the claim. g 



C Proof of Proposition [10 

(i) Set 

7 / \ - r . ( x + a \ (2x + a- l) 2 (x + l) 2 

\ x — 1 / x(x + a)(x — 1) 

Trivial but tedious computations show that the unique solution of h'(x) = such that 

x > 1 is x = x c which is a minimum, and that lim h{x) = lim h(x) = +00. Hence for 

x^l+ x— >+oc 

any > fc c , the equation V\(x, [x + a)/(x — l),x) = k has a unique solution in (l,x c ) and 
a unique solution in (x c , +00). Hence L^D C consists of two points, say p\ and p2- Recall 
that C is the curve of 2-periodic points of F, hence {pi,P2} is a 2-periodic orbit and (i) 
holds. 

To prove statement (ii) we will see that if k > k c the locus of non transversal intersections 
of Lfc with Q is the empty set. Consider the system G = and 

{Vi)x/{Vi) y = G x /Gy, (32) 
(ViW(Fi) 2 = G x /G z , (33) 

obtained by imposing {Wi||VG} Pi {G = 0}. 

The only positive solutions of equation (|32|) are given by the zeroes of mi := xaz + 
xaz 2 + 2ay 3 x + 3x 2 az + 5y 4 + 3y 5 + 2y 3 + 3x 2 az 2 + 2x 3 az + 2x 3 az 2 + xy 2 z — y 2 xz 2 + 6ay 2 z + 
4ay 3 z + 2xy 3 z + 2ay 2 x — 3y 2 x 2 z 2 + x 2 z + xz 2 + xz 3 + 3x 3 z — 2y 2 x 3 z 2 — 4y 2 x 3 + 2x 3 z 3 + 
2x 4 z 2 + 2x 4 z — 2y i x 2 — 2y 3 x 3 + 3x 2 z 3 + Ax 2 z 2 + 5x 3 z 2 — Ay 2 x 2 a — 2y 3 x 2 a — 2y 2 x 3 z — 2yx 2 a + 
2ayz — 2y 3 x 2 z — 2x 2 yz — 2yx 3 + 3y A x + 2y 3 z 2 + 5ay 4 + 3y 2 z 2 + 3y 2 z + 3ay 2 + 8y 3 z + 4y 3 x + 
8ay 3 + 3a 2 y 2 - Ay 3 x 2 - 3y 2 x 2 - yx 2 + a 2 y + 2a 2 y 3 + y 2 x + yz 2 + 5y 4 z - 7y 2 x 2 z, 
and its zeroes over are given by the the zeroes of n, where it satisfies mi = 91 G + n, 
and it is given by n := (2xz — 2x 3 z + 2xz 2 — 2x 3 z 2 ) y 2 + {Axz 2 + 2x 2 z + 2xaz 2 — 2x 3 z 
+Ax 2 z 2 + 2xz 3 + 2x 2 az 2 - 2x 3 z 2 - 2x A z 2 - 2x 4 z + 2x 2 z 3 + 2xaz + 2xz + 2x 2 az) y + 
2x 3 z 3 +2xaz+2xz 2 + 2xz 3 +4x 2 az 2 + 2x 3 az 2 +4x 2 az+2x 3 z 2 +4x 2 z 3 +2x 3 az+2xaz 2 +4x 2 z 2 . 

The only positive solutions of equation (|33|) are given by the zeroes of m2 := x — z and 
7713 := 2xaz + 2xaz 2 + 2x 2 az — y 3 + 2x 2 az 2 — ayx + 2xyz 2 — 2xy 2 z — ay 2 z — ay 2 x + 4x 2 z + 
4xz 2 + 3xz 3 + 3x 3 z + 2x 3 z 3 + 5x 2 z 3 + 9x 2 z 2 + 5x 3 z 2 — ya — ayz + 2x 2 yz + 2yz 2 x 2 — xy — 

q OO O OQ Q 00 O O O 

yz — y — y z — 2y z — ay — y z — y x — y x — yx — 2y x — yz + xz. 

Suppose that 7712 = so z = x and r\(x,y,x) = —2x 2 (x + l) 2 (l + x + y)(—a — x — y + xy), 
hence the positive solutions of n(x,y,x) = are given by y = (x + a)/(x — 1), thus the 
points in C But as mentioned above £P\Q = (x c , x c , x c ). 
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To study the zeroes of TO3 in Q, we consider the zeroes of r^, satisfying mi = q^G +ri- 
But it is given by r2 := 2xz{x + l)(z + l)(a + x + y + z + xz), and therefore = has not 
positive solutions. 

In summary, if k > k c , then VVi is never parallel to VG over Q and hence iti Q. 
(iii) Recall that Q is defined by the equation 

G = -y 3 - {x + z + a + l)y 2 - (x + z + a)y + xz{x + l)(z + 1) = 0. 

By applying Descartes' Rule on Q we obtain that for all x > and z > there exist a 
unique y(x,z) > solution of G = 0. Consider the function v(x, z) := Vi(x,y(x, z), z). 
Now the proof is done in two steps: (I) The only singular point of v is (x c ,x c ), which is a 
minimum. (II) Each level curve v(x, z) = k > k c is a closed curve surrounding (x c ,x c ). 
Step I: To find the singular points of h(x, z) we look for the solutions of system 



such that x > and z > 0. The only factors in v x and v z giving rise to such solutions are 

m := xaz + xaz 2 + 2ay 3 x + 3x 2 az + 5y 4 + 3y 5 + 2y 3 + 3x 2 az 2 + 2x 3 az + 2x 3 az 2 + xy 2 z — 
y 2 xz 2 + 6ay 2 z + 4ay 3 z + 2xy 3 z + 2ay 2 x - 3y 2 x 2 z 2 + x 2 z + xz 2 + xz 3 — 2y 3 x 2 a — Ay 2 x 2 a + 
3x 3 z + 2x 3 z 3 + 2x 4 z 2 + 2x 4 z + 3x 2 z 3 + 4x 2 z 2 + 5x 3 z 2 - 2y 2 x 3 z 2 - \y 2 x 3 + 3y 4 x + 2y 3 z 2 + 
5ay 4 + 3y 2 z 2 + 3y 2 z + 3ay 2 + 8y 3 z + 4y 3 x + 8ay 3 + 3a 2 y 2 — 4y 3 x 2 — 3y 2 x 2 — 2yx 3 — yx 2 + a 2 y + 
2a 2 y 3 + y 2 x + i/z 2 + 5y 4 z - 7y 2 x 2 z + 2ayz - 2yx 2 a - 2y 2 x 3 z - 2y 3 x 2 z - 2x 2 yz - 2y 4 x 2 - 2y 3 x 3 



n := xaz + 3xaz 2 + Aay 3 x + x 2 az + by 4 + 3y 5 + 2y 3 + 3x 2 az 2 + 2xaz 3 + 2x 2 az 3 + 2x 2 z 4 + 2xz 4 - 
2y 2 x 2 z 3 - 2y 2 xz 3 - 2yz 2 a - Ay 2 z 2 a - 2y 3 z 2 a - 2y 3 xz 2 + 2ayx - 2xyz 2 + xy 2 z - 7y 2 xz 2 + 
2ay 2 z + 2ay 3 z + 2xy 3 z + 6ay 2 x - 3y 2 x 2 z 2 + x 2 z + xz 2 + 3xz 3 + x 3 z + 2x 3 z 3 + 5x 2 z 3 + 4x 2 z 2 + 
3x 3 z 2 + 5y 4 x - 4y 3 z 2 + hay 4 - 3y 2 z 2 + y 2 z + 3ay 2 + 4y 3 z + 8y 3 x + 8ay 3 + 3a 2 y 2 + 2y 3 x 2 + 
3y 2 x 2 + yx 2 + a 2 y - 4y 2 z 3 - 2yz 3 - 2y 4 z 2 - 2y 3 z 3 + 2a 2 y 3 + 3y 2 x - yz 2 + 3y 4 z - y 2 x 2 z, 
respectively. Here y denotes y(x,z). 

The common zeroes of m and n in Q are given by the zeroes of the functions r and 
s respectively, where m = pG + r and n = qG + s for some polynomials p and q. These 
functions are 



The only positive solutions of r = and s = are given by (x, (x + a)/(x — l),x), which 
are the points of C. Since C n Q = (x c , x c , x c ), the proof of (I) is finished. 



< 



V X = (VI). + (Vl)y^ = 0, 

v z = (Vi) z + (Vl) y / = 0, 




and 



r = 2xz (x + 1) (1 + x + y) (z + 1) (a + y + z - xy), 
s = -2xz (x + 1) (1 + y + z) (z + 1) (-a - x — y + yz). 
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Step II: Note the following facts: from statement (ii), Q iti L^; the level curves v(x, z) = k, 
are defined by L/% n Q and thus for analytic equations and for k > k c they have no critical 
points (since otherwise the hamiltonian vector field —v z d x +v x d z would have another critical 
point than (x c ,x c ) in contradiction with Step I); The sets are compact. From all them 
we conclude that for a fixed k > k c , each level set of v(x,z) = k is diffeomorphic to a 
finite union of closed curves. To prove that indeed it is formed by an unique closed curve 
it suffices to show that the function x — > v(x,x c ) is monotonic in (x c ,+oo). To see this 
we will prove that it has the unique critical point x = x c . To this end note that the only 
positive solutions of r = are given by the factor a + y + z — xy, hence x = (a + y + z)/y. 
Now G{(a + y + x c )/y, y, x c ) = (y — x c ) qi{y)/y 2 , where is a degree four polynomial in y 
without positive solutions. Thus, as we wanted to SGGj X — X q IS the unique critical point of 
x — > v(x, x c ) and (iii) follows. gj 
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